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LIGHTNING  DATA  ANALYSIS 


SECTION  I.  INTRODUCTION 

'Phis  Final  Report  on  Contract  F3 36 1 5-8 1 -C-34 1 0  is  divided  in¬ 
to  three  parts:  Part  1  (SECTION  II)  discusses  our  analyses  of  cor¬ 
related  airborne  and  ground-based  electromagnetic  data  obtained 
during  the  AFWAL/FIESL  lightning  characterization  program  in  South 
Florida;  Part  2  (SECTION  III)  is  concerned  with  the  calculations 
of  lightning  return-stroke  electric  and  magnetic  fields  at  flight 
altitudes;  and  Part  3  (SECTION  IV)  includes  the  specification  of 
a  lightning  test  standard  and  a  discussion  of  the  validity  of 
deriving  lightning  currents  from  electric  and  magnetic  fields 
measured  remote  from  the  lightning.  This  part  also  discusses  some 
calculations  of  electric  field  frequency  spectra  and  a  comparison 
of  lightning  and  an  NEMP  produced  by  an  exoatmospher ic  burst. 

SECTION  II.  CORRELATED  AIRBORNE  AND  GROUND-BASED  LIGHTNING  ELECTRO¬ 
MAGNETIC  DATA 

Durinq  1979,  1980,  and  1981  the  AFWAL/FIESL  directed  a  program 
designed  to  characterize  airborne  lightning  electric  and  magnetic 
fields  (see  "Airborne  Lightning  Characterization”,  AFWAL-TR-83-30 1 3 , 
January  1983,  by  P.  L.  Rustan,  B.  P.  Kuhlman,  A.  Serrano,  J.  Reaser, 
and  M.  Risley) .  A  WC-130  aircraft,  instrumented  for  electric  and 
magnetic  field  measurements,  was  flown  in  South  Florida  in  the 
vicinity  of  a  network  of  ground-based  stations  that  provided  wide¬ 
band  electric  fields  at  ground  level  and  data  that  could  be  used  to 
determine  the  location  of  lightning  VHF  sources.  A  thorough  analysis 
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of  the  voluminous  data  obtained  will  take  many  years.  We,  in 
conjunction  with  AFWAL  personnel  ,  have  completed  a  survey  analysis 
of  ten  lightning  events  from  the  South  Florida  study  and  a  detailed 
analysis  of  one  event  that  occurred  at  17:09:40  EDT  on  July  16, 

1981.  Appendix  A  contains  a  scientific  paper  that  describes  the 
results  of  our  analysis  of  the  17:09:40  event.  This  flash  is 
typical  of  most  of  the  lightning  events  in  that  the  results  obtained 
and  new  questions  raised  are  similar  to  those  for  other  events. 

For  the  17:09:40  event,  a  two-stroke  flash  with  separate  chan¬ 
nels  to  ground,  lightning  channel  reconstructions  were  possible  for 
both  channels  to  ground  via  the  VHF  time-of -arrival  system.  Elec¬ 
tric  and  magnetic  fields  were  recorded  on  the  WC-130  and  electric 
fields  were  recorded  at  the  qround .  The  following  conclusions  can 
be  drawn  from  the  detailed  analysis  of  the  17:09:40  flash  and  sur¬ 
vey  analysis  of  other  events:  Maximum  rates-of -change  of  airborne 
electric  fields  from  the  last  few  stepped  leader  pulses  and  from  the 
fast  field  transitions  of  return  strokes  are  about  3  V/m  ys  normal¬ 
ized  to  100  km.  These  values  are  an  order  of  magnitude  smaller  than 
those  reported  by  Weidman  and  Krider  (1980)  and  Weidman  (1982)  for 

leader  and  return  stroke  pulses  measured  at  ground  level  over  salt 

dE 

water.  Figure  1  shows  a  histogram  of  maximum  /dt  values  measured 
at  ground  level  over  salt  water  and  normalized  to  100  km  for  the  fast 
field  transitions  of  return  strokes.  The  mean  is  33  V/m  ys.  It  is 
crucial  to  understand  why  the  airborne  electric  field  rates-of-chanqe 
are  so  slow.  One  possible  explanation  that  needs  further  examination 
is  that  since  the  stepped  leader  and  return  stroke  fields  are  gener¬ 
ated  near  the  ground,  the  electromagnetic  wave  may  be  coupled  to  or 
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Figure  1.  A  Histogram  of  the  Measured  Maximum  AE/At  Values 
for  Return  Strokes  Normalized  to  a  Range  of  100  km 
(from  Weidman,  1982). 
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reflected  from  the  non-perfectly-conductina  earth  and  thereby  produce 
a  depredation  of  the  airborne  risetime.  A  theoretical  analysis  is 
possible.  Additionally,  the  fields  at  the  aircraft  altitude  or  above 
should  be  examined  to  see  if  there  are  faster  rates-of-channe  that 
are  characteristic  of  pulses  produced  far  above  ground.  A  search  of 
the  records  for  these  kinds  of  data  is  outside  the  scope  of  the  pre¬ 
sent  study. 

Clectric-field  risetimes  measured  on  the  aft  uoper  fuselaae  of 
the  WC-130  are  always  faster,  up  to  0.5  ;;s  faster,  than  those  mea¬ 
sured  on  the  forward  upper  fuselaae.  Further,  the  field  amplitudes 
on  the  aft  upper  fuselaae  are  a  factor  of  3  to  4  smaller  than  on  the 
forward  upper  fuselage.  The  model  study  by  Electro  Maanetic  Applic¬ 
ations,  Inc.  (EMA)  of  the  electromagnetic  response  of  the  WC-130 
aircraft  to  the  e lectromaanctic  signals  originatinn  at  the  known 
locations  of  the  stepped  leaders  and  return  strokes  of  the  17:09:40 
event  show  the  risetimes  of  the  two  sensors  to  be  identical.  More¬ 
over,  each  antenna  response  is  essentially  the  incident  field  times 
a  scale  factor.  A  possible  explanation  for  the  observed  difference 
in  the  two  sensors  is  that  on  the  actual  aircraft  there  was  a  thin- 
wire  antenna  above  the  aft  upper  fuselage  antenna  which  both  could 
have  shielded  it  and  could  have  caused  resonances  that  affected  the 
measured  risetime.  It  was  not  reported  until  very  recently  that 
there  was  such  a  thin-wire  antenna  on  the  WC-130,  so  the  antenna  was 
not  modeled  in  the  EMA  study.  Future  work  should  include  modeling  of 
the  exact  configuration  of  the  aircraft  including  all  wire  antennas. 
The  measured  fields  on  the  aircraft  show  a  pronounced  resonance  at 
about  3.7  MHz  which  is  not  produced  by  the  model  illumination  of  the 
WC-130  with  a  spherical  electromagnetic  field.  The  model  calculations 
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show  a  relatively  small  resonance  at  about  that  frequency.  In  order 
to  understand  better  the  generation  of  resonances,  the  model  calcula¬ 
tions  should  be  expanded  to  include  field  components  in  the  direction 
of  propagation  which  may  be  more  efficient  excitors  of  resonances. 

The  effects  of  angle  of  incidents  of  the  incoming  el ectromacmet  ic  wave 
on  the  aircraft  response  should  also  be  studied. 

SECTION  III  AIRBORNE  ELECTRIC  AMD  MAGNETIC  FIELD  CALCULATIONS 

The  results  of  our  calculations  of  lightning  return  stroke  fields 
above  ground  are  given  in  the  scientific  paper  in  Appendix  B.  Two 
sets  of  calculations  are  presented:  those  for  which  the  initial  pea'--, 
current  propagates  up  the  channel  unattenuated  and  those  for  which,  the 
peak  decreases  exponentially  with  height  with  a  decay  length  of 
km.  Recent  work  by  Jordan  and  Uman  (1983)  has  shown  that,  while  the 
return  stroke  light  decays  exponentially  with  height  with  about  a  1.0 
km  decay  length,  the  relation  between  light  and  current  is  such  that, 
in  all  probability,  the  current  decays  much  more  slowly  with  height 
than  does  the  light.  Thus,  the  fields  to  be  expected  above  ground 
should  be  between  the  two  cases  presented  in  Appendix  B. 

The  results  in  Appendix  B  are  primarily  concerned  with  overall 
field  waveshapes  and  no  data  are  given  on  nsetimes.  The  risctimes  in 
the  air  above  a  perfectly  conductina  ground  are  determined  by  both,  the 
risetime  of  the  line-of -sight  wave  and  the  delayed  risetine  of  the 
reflected  wave  from  the  ground.  For  each  small  channel  section  that 
radiates  a  field  the  propagation  delay  is  different,  and  thus  the 
total  risetime  observed  above  a  perfectly-conductina  ground  can  only 
be  determined  by  a  relatively  complex  calculation.  In  view  of  the 
fact  that  the  risetimes  measured  on  the  WC-130  aircraft  are 


substantially  slower  than  ground-based  measurements  made  over  salt 
water,  as  noted  in  the  previous  section,  the  calculation  of  airborne 
risetimes,  oven  for  the  case  of  a  perfectly  coniuctinq  earth,  would 
appear  to  be  desirable. 

SECTION  IV.  LIGHTNING  TEST  STANDARD 

Tables  IA  and  IIA  of  the  scient  ic  paper  in  Appendix  C  give 
our  recommendation  for  a  lightning  1  urn  stroke  test  standard  for 
average  first  and  subsequent  return  ces,  as  well  as  a  comparison 

of  lightning  and  NEMP .  We  recommend,  as  noted  in  Appendix  III,  that 
severe  strokes  be  modeled  by  scaling  the  current  amplitudes  for  aver¬ 
age  strokes  up  by  a  factor  of  5.  We  have  chosen  to  use  current  wave¬ 
forms  derived  from  remote  field  measurements  rather  than  those  found 
from  direct  tower  measurement  (given  in  Tables  IB  and  IIB)  for  the 
reasons  discussed  in  Appendix  C  and  in  the  following.  The  best  avail¬ 
able  direct  current  measurements,  as  discussed  in  Appendix  C,  come 
from  Berger  and  Garbagnati  and  are  based  on  strikes  to  towers  on  two 
mountains  near  the  Swiss-Italian  border.  In  these  data,  the  risetimes 
of  first  return  strokes  are  considerably  slower  than  those  of  subse¬ 
quent  strokes  and  a  peak  current  derivative,  °'t/dt,  of  1  x  10:1  A/s 
occurs  in  about  1?,  of  the  subsequent  strokes.  However ,  in  a  South 
African  study,  as  noted  in  Appendix  C,  a  L,:i/dt  of  1.8  x  10  !l  A/s  was 
measured  for  one  lightning  strike  to  a  tower  on  relatively  flat  ground 
in  a  small  sample  of  flashes.  This  case  is  apparently  the  larcrest 
°^^dt  that  has  been  measured  directly.  Later  measurements  on  the  same 
South  African  tower  (Erikkson,  1982)  using  a  waveform,  digitizer  showed 
4  out  of  5  subsequent  strokes  in  a  six  stroke  flash  (out  of  a  sample 
of  3  multiple-stroke  flashes)  had  risetimes  less  than  0.2  ;.s,  the 
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sample  t  ia.e ,  and  the  ^  /dt  ~  qrouter  than  0.08,  1.5,  i.O,  am: 

times  1 0  1  ’  A/s.  Whether  currents  measured  on  tower.-,  are  •  r  u  1 
resentative  of  the  currents  in  the  1 ightning  channel  above  n 
or  of  the  currents  that  would  1  low  throuqh  an  aircrait  aiio'.v 
is  not  known.  The  shape  of  the  tower  current,  particular. y  t 
the  first  return  stroke  in  a  flash,  is  not  consistent  with 
trie  and  magnetic  fields  produced  by  normal  1 iahtninq  to  nrcu. 
(Weidman  and  Krider,  1978).  Unfortunately,  there  are  no  situ 
measurements  of  the  electromagnetic  fields  and  currents,  dar  i  r. 
quent  strokes  in  rocket- triggered  f  lashes ,  am;  t  hey  sv«  ■  :a  • 
measurements  to  compute  return  stroke  ve  1  oei  t  ies  a.;;,  ;  t:;e  *.  •: 
aiven  in  Appendix  C.  Using  this  method,  the  Krone;.  .  btu  ;a<. 
velocity  of  1.3  x  10'  m's  with,  a  standard  C1GV  ldt  iO!l  v  ' !  v»  .  A 
using  magnetic  fields,  and  a  mean  of  1.7  x  10  rr/s  i  .1: 

deviation  of  0.43  x  10"  m/s  using  electric  fields  (u  ieux  oi  a 
Dejebari  et  al .  ,  1981).  These  velocity  c2e  t  o  rm  i  na  t .  on  s  a  re  c 
with  the  photographic  measurements  of  I  clone  and  nrvi  I  le  (  j  •< 
report  a  mean  of  0.96  x  10  m/s  for  first  strokes  within  1  hr 
ground  and  a  1.2  x  10  m/s  for  subsequent  strokes.  "hero:  or  ■ 
regard  the  French  measurements  on  triune  rod  !  ;  u:.t  u ;  u.u  as  .  i  v 
some  support  for  the  theory  aiven  in  Appendix  U.  It  is  inter 
to  note  that  the  10-90  1  risetime  of  the  Kroner.  Ci.nvnt  x 
shown  as  an  example  is  about  0.1  us  and  that  the  peak  current 
10  kA  (Fieux  et  al.,  1978)  ,  yield.inq  a  ^  1  dt  of  1  x  10  A/s. 

Return  stroke  currents  that  are  derived  from  n ensured  ft 
a  mean  maximum  c^/dt  of  about  1.5  x  10‘  1  A/s,  and  the  maximum, 
value-'  is  about  4  x  1 0 1  1  A/s  in  about  100  measurements.  There 
moan  maximum  <,l/cit  derived  from  fields  is  cqui  va  ie-r.t  t;  •re 
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in  the  tower  data. 


In  the  paper  in  Appendix  C,  we  have  assumed  that 
a  typical  lightning  has  maximum  uI/dt  of  1.5  x  lu n  A/s  and  peak  cur¬ 
rent  of  35  kA,  and  that  a  severe  lightning  has  a  maximum  ^*/dt  and 
peak  current  that  are  five  times  larger  than  the  eypical  liqhtning, 
i.e.  7.5  x  1011  A/s  and  175  kA,  respectively.  These  choices  for  a 
severe  lightning  have  been  criticized  because  they  associate  the 
largest  peak  current  with  the  greatest  C*I/dt.  We  shall  explore  the 
validity  of  these  choices  below. 

Figure  1  shows  the  submicrosecond  structure  of  a  typical  return 
.arokc  radiation  field  and  identifies  the  portion  just  prior  to  the 
peak  that  has  the  largest  uE/dt.  In  our  model,  ^‘"/dt  is  directly 
proportional  to  ^/dt  and  the  constant  of  proportionality  contains 
the  return  stroke  velocity  near  ground,  as  noted  in  Appendix  C.  A 
h  i  s  toe  ram  of  measured  ^f'/dt  values  normalized  to  100  km  for  the  fast 
field  transition  are  plotted  in  Figure  1  for  lightning  at  a  number  of 
distances  over  salt  water  (Weidman  and  Krider,  1980;  Weidman,  1982). 
i'hese  measurements  were  made  over  salt  water,  and  evidently  the  pro¬ 
pagation  distance  does  not  affect  the  measured  values.  The  mean  max- 
i  i  H 

imum  ‘/dt  during  the  fast  transition  is  33  V/m  ys  normalized  to  100 
km,  and  the  mean  10  to  90t>  field  risetime  is  90  nsec  during  the  fast 
transition. 

Figure  2  shows  the  relationship  between  the  maximum  l"/dt  ami 

d  F 

the  corresponding  AE  during  the  fast  transition.  The  values  of  '/dt 
and  E  do  appear  to  be  correlated,  and  this  implies  that  a  larae  cur¬ 
rent  peak  will  produce  a  large  <^I/dt  as  we  have  assumed  for  our 
severe  lightning.  On  the  other  hand,  only  8  out  of  108  points  in 
Figure  2  are  above  50  V/m  us,  and  these  have  a  larger  variation  in 
!■;  than  the  points  below  50  V/m  ys.  Therefore,  it  might  be  argued 


that  there  is  not  enouqh  ^E/dt  data  to  draw  a  firm  conclusion  about 
the  distribution  at  high  values  of  ^'"/dt.  It  has  been  suggested  that 
the  data  may  be  approaching  a  limit  at  about  75  V/m  us,  but  this  does 
not  appear  to  be  valid  in  view  of  the  small  number  of  measurements. 

It  has  also  been  suggested  that  the  data  above  50  V/m  ;is  may  be  pro¬ 
duced  by  a  different  process  than  the  data  below  this  value,  e.g.  by 
two  channels  radiating  simultaneously,  but  there  is  still  no  direct 
evidence  that  this  suggestion  is  valid. 

Figure  3  summarizes  all  the  available  data  on  the  values  of  the 

maximum  return  stroke  ^I/dt.  The  data  for  the  field-derived  ^/dt 

are  plotted  assuming  a  return  stroke  velocity  of  1  x  1 0 m/s.  The 

plotted  lines  show  where  these  data  would  fall  if  the  velocity  were 

either  1.4  x  10  or  0.6  x  10fl  m/s.  It  is  clear  from  this  figure  that 

d  I  1  s 

the  average  field-derived  /dt  correspond  to  the  maximum  values  of 
the  tower  measurements  for  subsequent  strokes  and  that  the  tower 
values  for  first  strokes  are  significantly  lower  than  those  for  sub¬ 
sequent  strokes. 

As  noted  earlier,  the  validity  of  our  model  relating  fields  and 
currents  is  supported  by  the  French  measurements  on  triggered  light- 
nina,  and  we  think  this  model,  which  assumes  that  an  upward  propaga¬ 
ting  current  pulse  is  associated  with  the  return  stroke  wavefront, 
is  the  best  that  is  currently  available.  An  alternate  model,  which 
assumes  that  a  spatially  uniform  but  time-varying  current  propagates 
upward,  the  so-called  Bruce-Golde  model,  yields  a  field-derived 
^/dt  that  is  within  a  factor  of  2  of  that  found  with  our  model. 
Essentially  all  of  the  high  frequency  content  of  the  field  is  deter¬ 
mined  by  the  current  rise  to  peak  and  the  current  fall  just  after 
peak,  so  the  validity  of  the  return-stroke  current  model  after  the 


Figure  3.  The  Distribution  of  the  Maximum  /dt  in  Return  Strokes  Derived 
from  Field  Measurements  {from  Weidman,  3902). 


peak  is  oi  secondary  importance. 


It  has  also  been  suggested  (Liman  et  al  .  ,  1973  ;  Weidman  and 
Krider,  1978)  that  the  initial  first  stroke  field  may  be  produced 
by  currents  propagating  both  upward  and  downward  from  the  junction 
between  the  upward  and  downward  leaders.  This  effect  would  lower 
our  field-derived  dI/dt  by  a  factor  of  2;  but  such  an  effect  should 
not  occur  in  subsequent  strokes  and  these  are  observed  to  have  about 
the  same  dE/dt  and  hence  dI/dt  as  first  strokes. 

As  an  extension  of  the  calculations  in  Appendix  C,  Figure  4 
shows  electric  field  spectra  for  an  average  first  stroke  at  dis¬ 
tances  between  50  m  and  10  km.  The  dashed  lines  in  Figure  4  show 
the  spectra  of  just  the  radiation  field  term  so  the  contribution 
of  the  electrostatic  and  induction  fields  can  be  evaluated. 

An  aircraft  in  flight  probably  will  not  encounter  the  return 
strokes  current  characteristic  of  ground  level  considered  above; 

dp 

but,  on  the  other  hand,  the  maximum  "/dt  values  in  cloud  pulses 
and  leader  steps  and  the  associated  amplitude  spectra  above  10'  Hz 
are  very  similar  to  those  of  return  strokes  (Weidman  et  al.,  1981). 
Although  we  do  not  yet  have  a  model  for  these  processes  in  which 
we  are  confident,  the  available  measurements  imply  that  the  maximum 
current  derivatives  in  these  processes  are  comparable  to  return 
strokes.  Therefore,  we  expect  that  the  hazards  from  the  high  fre¬ 
quency  components  of  cloud-discharges  may  well  be  similar  to  return 
strokes  near  the  ground. 

Since  the  actual  lightning  channel  is  tortuous,  it  might  be 
expected  that  the  frequency  spectra  shown  in  Figure  4  miqht  be  af¬ 
fected  by  that  tortuosity.  We  argue  now  that  this  is  not  the  case 
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at  close  range  where  the  large  field  magnitudes  could  adversely 
affect  an  aircraft.  Levine  and  Meneghini  (1978a)  have  used  a  simple 
current  model  to  calculate  the  fields  which  are  radiated  by  a  tortu¬ 
ous  channel  and  have  shown  that  the  tortuosity  increases  the  "jag¬ 
gedness"  of  a  time-domain  waveform  and  increases  the  spectral  ampli¬ 
tude  above  10 5  Hz  by  about  20  db.  We  have  repeated  their  calcula¬ 
tions  for  both  distant  and  close  (50  m)  fields  for  a  first  stroke 
that  has  the  current  parameters  given  in  Appendix  C,  Table  1A.  The 
channel  tortuosity  is  that  given  in  Figure  2  of  Levine  and  Meneghini 
(1978b).  The  results  are  shown  in  Figure  5,  and  it  is  clear  that 
tortuosity  does  not  appreciably  alter  the  spectrum  of  the  electro¬ 
static  or  induction  components  which  dominate  the  fields  at  close 
ranges.  The  change  in  the  radiation  field  spectrum  with  tortuosity 
is  an  increase  of  about  10  db  above  10s  Hz.  These  calculations  are 
critically  dependent  on  the  channel  current  waveform  and  the  assumed 
tortuosity.  On  the  other  hand,  the  time-domain  waveforms  for  the 
simulated  tortuosity  are  much  more  "jaqged"  that  the  experimental 
data,  which  for  subsequent  strokes  are  actually  quite  smooth,  so  the 
effects  of  tortuosity  may  not  be  nearly  as  large  as  these  calculations 
would  indicate.  In  fact,  most  of  the  frequency  content  above  10 6  Hz 
in  the  measured  time-domain  fields  from  first  and  subsequent  strokes 
is  produced  within  1  usee  or  so  of  the  peak  field;  and  this  implies 
that  most  high  frequencies  are  radiated  at  a  time  when  the  stroke  is 
within  a  few  hundred  meters  of  ground  and  prior  to  the  time  when 
tortuosity  can  play  a  significant  role.  Why  subsequent  stroke 
field  waveforms  are  smooth  when  photographed  channels  appear  to  have 
considerable  tortuosity  is  not  clear.  Currently,  there  are  studies 
under  way  at  the  University  of  Arizona  to  measure  tortuosity  and 


branching  in  real  channels,  and  in  the  future  these  will  be  coupled 
with  calculations  of  fields  by  the  University  of  Florida.  We  hope 
that  these  future  studies  will  resolve  this  apparent  contradiction. 
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ABSTRACT 


We  have  reduced  and  analyzed  the  data  from  a  two-stroke  lightning 
flash  to  ground  which  occurred  in  South  Florida  on  July  16,  1981,  within 
a  network  of  four  ground  stations  instrumented  for  VHF  measurements, 
about  15  km  from  a  ground  station  instrumented  for  wide-band  electric 
field  measurements,  and  within  10  km  of  a  WC-130  aircraft  operating 
at  5.2  km  and  instrumented  for  wide-band  electric  and  magentic  field 
measurements.  The  four— station  ground-based  VHF  measurements  allow 
a  reconstruct!  <n  of  the  geometry  of  the  flash,  which  was  composed  of 
two  separate  channels  to  ground.  Electric  field  system  bandwidth  for 
the  ground  measurement  was  from  0.02  Hz  to  about  2  MHz;  electric  and 
magnetic  field  system  bandwldths  on  the  aircraft  extended  to  20  MHz. 
Ground— baaed  and  airborne  measurements  of  fields  are  presented  and 
shown  to  be  consistent  with  one  another. 
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DURING  1979,  I960,  AND  1981  THE  AIR  FORCE 
Wright  Aeronautical  Laboratories  directed  a 
prograa  designed  to  characterise  airborne 
lightning  electric  anu  magnetic  fields. 

A  WC-130  aircraft  instrumented  for  electric 
and  magnetic  field  measurements  was  flown  in 
South  Florida  in  the  vicinity  of  a  network 
of  ground-based  stations  which  provided 
electric  field  at  ground  level  and  data  from 
which  the  location  of  lightning  VHF  sources 
could  be  determined.  Extensive  data  were 
obtained.  These  will  take  many  years  to 
analyze  fully.  In  this  paper  we  briefly 
describe  the  instrumentation  used  during  the 
1981  measurement  season  and  Illustrate  the 
potential  of  the  data  base  by  presenting  an 
analysis  of  data  from  one  two-stroke  light¬ 
ning  flash  to  ground  which  occurred  at 
17:09:40  EDT  on  July  16,  1981. 

AIRBORNE  MEASUREMENT  SYSTEMS 

The  WC-130  aircraft  is  about  30  m  from 
nose  to  tall  and  about  41  m  in  wingspan. 
Aircraft  resonances  are  expected  at  half  and 
Integer  multiples  of  9.9  and  7.4  MHz.  The 
airborne  Instrumentation  had  an  upper  fre¬ 
quency  response  limit  of  about  20  MHz  so  that 
some  of  these  resonance  effects  could  be 
observed.  Three  basic  types  of  sensors, 
described  in  Baum  et  al.  (l)*twere  used: 

(a)  plates  to  measure  the  component  of  the 
electric  field  intensity  perpendicular  to 
aircraft  surfaces,  (b)  loops  to  measure  the 
magnetic  field  intensity  parallel  to  aircraft 
surfaces,  and  (c)  loops  to  measure  current 
densities  flowing  in  aircraft  surfaces  by 
sensing  the  magnetic  field  associated  with 
those  current  densities.  The  latter 
two  sensors  have  essentially  similar  principles 
of  operation.  A  total  of  eleven  sensors  were 
used  on  the  WC-130  in  1981.  Electric  field  was 
measured  on  the  forward  upper  fuselage,  aft 
upper  fuselage,  aft  lower  fuselage,  and  left 
wing  tip.  Both  horizontal  components  of  the 
magnetic  field  were  measured  on  the  forward 
upper  fuselage.  Skin  current  density  was  mea¬ 
sured  on  the  top  and  bottom  of  each  wing  and 
on  the  aft  upper  fuselage.  All  measured  quan¬ 
tities  were  continuously  recorded  on  instrumen¬ 
tation  tape  with  an  upper  frequency  response 
limit  of  about  2  MHz.  In  addition,  the  deriva¬ 


tives  of  the  measured  quantities  were  sampled 
at  20  ns  intervals  for  time  blocks  of  160  ps. 

Such  blocks  of  data,  with  an  effective  upper 
frequency  response  limit  of  about  20  MHz,  were 
acquired  at  a  rate  of  twice  a  second,  the  data 
block  being  initiated  in  a  pre-trigger  mode  by 
an  incoming  signal  exceeding  a  pre-set  thresh- 
hold. 

GROUND-BASED  MEASUREMENT  SYSTEMS 

Ground-based  electric  field  measurements  are 
essential  to  proper  interpretation  of  the  air¬ 
borne  data  since  considerable  information  exists 
on  the  characteristics  of  the  fields  observed 
at  ground  level  and  the  relation  of  those 
fields  to  their  sources,  whereas  such  informa¬ 
tion  is  not  available  for  airborne  fields.  The 
ground-based  electric  field  system  was  similar 
to  that  described  in  Beasley  et  al.  (2).  The 
fields  were  recorded  on  eight  channels  of  an 
instrumentation  tape  recorder  with  a  bandwidth 
in  the  FM  mode  of  0.02  Hz  to  500  kHz  and  in  the 
direct  mode  of  400  Hz  to  2  MHz.  A  variety  of 
gains  allowed  the  measurement  of  fields  between 
4  V/m  and  40,000  V/m.  Fig.  1  shows  the  overall 
experimental  setup  Including  the  location  of 
the  trailer  that  housed  the  electric  field 
system. 

The  VHF  source  location  system  comprised 
four  VHF  stations  located  about  20  km  apart  as 
shown  in  Fig.  1.  The  VHF  radiation  at  each 
station  was  (a)  detected  with  an  omnidirectional 
antenna,  (b)  passed  through  a  filter  with  a 
center  frequency  of  63  MHz  and  a  bandwidth  of 
6  MHz,  (c)  log  amplified,  (d)  envelope  detected, 
and  (e)  recorded  on  a  modified  version  of  the 
RCA  VCT  201  Video  Cassette  Recorder.  The  sys¬ 
tem  allows  VHF  locations  from  the  measurement 
of  the  difference  in  the  time  of  arrival  of  a 
given  pulse  at  the  four  stations  as  explained 
in  Rustan  et  al.  (3)  and  Proctor  (4).  The  time 
correlation  necessary  for  this  measurement,  about 
0.1  us,  was  accomplished  by  using  WWV  for  crude 
time  correlation  and  the  vertical  and  horizontal 
sync  pulses  from  WINK-TV  in  Fort  Myers  (shown 
in  Fig.  1)  for  fine  time  correlation.  The  hor¬ 
izontal  sync  pulses  have  a  rate  of  one  each  63 
us. 


♦Numbers  in  parentheses  designate  References 
at  end  of  paper. 
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A  conceptual  sketch  of  the  lightning 
channels  of  a  two-stroke  flash  occurring  on 
July  16,  1981  at  17:09:40  EDT  is  shown  In  Fig. 

1.  The  sketch  is  based  on  the  VHF  time  of 
occurrence  and  location  of  the  VHF  radiation 
sources  shown  in  plan  view  in  Fig.  2a  and  look¬ 
ing  north  in  Fig.  2b.  The  location  and  orien¬ 
tation  of  the  WC-130,  which  was  flying  at  5.2 
km,  is  also  shown  in  both  figures. 

Both  strokes  appeared  to  originate  from 
about  the  same  region,  but  the  second  went  to 
a  different  ground  strike  point,  about  5  km 
north-west  about  250  ms  after  the  first  stroke. 

The  first  VHF  radiation  sources  start 
about  50  ms  before  the  first  return  stroke,  at 
an  altitude  of  about  7  km,  that  is,  2  km  above 
the  level  of  the  WC-130,  and  about  7  km  east 
and  3  km  south  of  it.  The  source  locations  then 
spread  up  and  down  about  1  km  and  east  to  about 
8  km  in  about  5  ms.  During  the  last  5  ms  before 
the  return  stroke,  the  source  locations  are  at 
an  altitude  between  4  km  and  1  km,  from  7  to 
10  km  east  and  from  2  to  6  km  south  of  the 
WC-130.  The  ground  strike  point  of  the  first 
stroke  appears  to  be  between  7  and  8  km  east 
and  2  to  3  km  south  of  the  WC-130.  For  about 
0.5  ms  after  the  return  stroke,  VHF  sources 
appear  between  3  and  8  km  altitude,  7  to  10  km 
east  and  1.5  to  8  km  south  of  the  WC-130. 

About  200  ms  later,  VHF  sources  become 
active  for  0.5  ms  between  6  and  3  km  altitudes, 

5  to  9  km  east  and  2  to  3  km  south  of  the 
WC-130.  Then,  30  ms  later,  for  about  1  ms,  VHF 
sources  appear  from  5  km  down  to  1.5  km  alti¬ 
tudes,  4  to  5  km  east  and  2  to  5  km  south  of 
the  WC-130.  The  strike  point  appears  to  be 
about  5  km  east  of  the  WC-130. 

Figs.  3,  4,  and  5  show  the  vertical  elec¬ 
tric  field  at  the  ground  station,  the  airborne 
vertical  electric  field  on  the  aft  upper  fuse¬ 
lage  (AUF  in  Fig.  1),  and  the  magnetic  field  in 
the  direction  of  the  fuselage  as  measured  on  the 
forward  upper  fuselage  (FUF  on  Fig.  1),  respec¬ 
tively,  for  the  first  stroke  in  the  flash. 

The  stroke  which  produced  the  electric 
field  in  Fig.  3  was  at  a  range  of  16  km  from 
the  ground  station  and  8  km  from  the  aircraft. 
The  airborne  field  magnitudes  in  Figs.  4  and  5 
are  not  corrected  for  field  distortion  by  the 
aircraft.  The  stepped  leader  pulses  which  pre¬ 
cede  the  return-stroke  transition  and  the  first 
ten  microseconds  or  so  of  the  return-stroke 
field  are  essentially  radiation  field  at  these 
ranges;  and  the  fields  on  and  above  the  ground 
are  expected  to  have  essentially  the  same  shape 

(5) ,  as  the  results  in  Figs.  3,  4,  and  5  con¬ 
firm.  After  about  10  us,  the  return-stroke 
electric  fields  show  an  electrostatic  component 
which  the  magnetic  field  does  not  possess  (5) 

(6) .  Additionally,  the  low-frequency  cut  off 
of  the  system  used  to  obtain  the  magnetic  field 
shown  in  the  figure  decays  slightly  faster  than 
the  actual  field. 


The  first-stroke  electric-field  intensity 
measured  at  the  ground  station  has  an  Initial 
peak  value  of  50  V/m,  or  8.0  V/m  normalized  to 
100  km,  a  peak  field  typical  of  Florida  light¬ 
ning  (7).  The  comparable  fiel  values  at  the 
WC-130  are  110  V/m,  or  a  normalized  8.8  V/m, 
on  the  forward  upper  fuselage,  and  32  V/m,  or 
a  normalized  2.6  V/m,  on  the  aft  upper  fuselage. 
The  second  stroke  peak  field  measured  on  the 
ground  was  14  V/m,  and  the  stroke  was  at  a 
range  of  14  km,  resulting  in  a  normalized  field 
of  2.0  V/m,  a  relatively  small  value  for  Florida 
return  strokes  (7).  The  airborne  second  stroke 
fields  were  comparably  small  and  difficult  to 
make  any  measurements  on  other  than  amplitude. 

In  the  first  stroke  field  records  the  ratios  c 
the  stepped-leader  pulse  heights  to  the  return 
stroke  peak  are  essentially  the  same,  on  aver ve 
about  0.1,  and,  as  expected,  the  stepped- leaser 
pulses  occur  at  the  same  times  before  the  re¬ 
turn  stroke  on  all  three  records,  as  can  be 
seen  in  Figs.  3,  4,  and  5.  The  zero-to-peak 
rise-time  of  the  first  return  stroke  measured 
at  the  ground  is  about  3.0  us,  also  consistent 
with  typical  values  measured  in  Florida  (7) , 
with  airborne  values  of  2.9  us  at  the  forwarc 
upper  fuselage  and  2.6  ys  at  the  aft  upper 
fuselage.  Stepped  leader  pulses  have  zero-to- 
peak  rise-times  on  the  ground  of  about  1  ps 
and  full-widths  at  the  pulse  base  of  about  -  _  . 
Comparable  airborne  values  are  0.8  ps  and  ... 
ps  at  the  forward  upper  fuselage  and  0. 3  jk  a no 
1.0  ps  at  the  aft  upper  fuselage.  All  measured 
rise-times  are  well  within  system  limits.  The 
rise-times  on  the  ground  are  expected  to  be 
longer  than  in  the  air  because  of  the  effects 
of  groundwave  propagation  involving  a  non- 
perfectly-conducting  earth,  as  discussed  In 
Lin  et  al.  (7),  Uman  et  al.  (8),  and  Weidman 
and  Krider  (9).  The  reason  that  the  rise-tirr.es 
at  the  aft  upper  fuselage  are  faster  than  then  .■ 
at  the  forward  upper  fuselage  is  not  known,  :  .. 
aircraft  resonances  may  contribute  to  this 
effect.  Wing  and  fuselage  resonances  are 
excited  by  the  airborne  horizontal  electric 
field  which  is  the  dominant  field  within  about 
1  km  of  a  return  stroke  and  is  comparable  to 
the  vertical  field  near  10  km  (5).  The  leader 
pulse  rise-times  are  somewhat  slower  than  the 
typical  values  for  10  to  90  percent  of  0.1  ..s 
reported  in  Weidman  and  Krider  (10)  for  light¬ 
ning  over  salt  water. 

Maximum  rates-of  change  of  airborne  elec¬ 
tric  field  for  both  leaders  pulses  and  return 
strokes  were  the  same,  about  40  V/r.  ps,  or 
3.2  V/m  us  normalized  to  100  km.  These  are  to 
be  compared  with  the  normalized  mean  of  30  V'rr. 
us  for  return  strokes  and  21  V/m  „s  for  leader 
pulses  reported  in  Weidman  and  Krider  (9)  1,10) 
for  lightning  over  salt  water. 


51-3 

21 


ACKNOWLEDGEMENT 


This  research  was  supported,  in  part,  by 
the  Wright  Aeronautical  Laboratories,  Wright 
Patterson  Air  Force  Base,  Ohio  under  contract 
F33615-79-K-0177. 

REFERENCES 

1.  C.E.  Baum,  E.L.  Breen,  J.P.  O'Neill, 
C.B.  Moore,  and  G.D.  Sower,  "Electromagentic 
sensors  for  general  lightning  application," 
in  Lightning  Technology,  NASA  Conference 
Publication  2128,  FAA-RD-80-30,  1980. 

2.  W.H.  Beasley,  M.A.  Uman,  and  P.L. 
Rustan,  Jr. ,  "Electric  fields  preceding  cloud- 
to-ground  lightning  flashes,"  J.  Geophy s.  Res,, 
87,  4883-4902,  1982. 

3.  P.L.  Rustan,  Jr.,  M.A.  Uman,  D.G. 
Childers,  W.H.  Beasley,  and  C.L.  Lennon, 
"Lightning  source  locations  from  VHP  radiation 
data  for  a  flash  at  Kennedy  Space  Center," 

J.  Geophys.  Res.,  85,  4893-4903,  1980. 

4.  D.E.  Proctor,  "A  hyperbolic  system  for 
obtaining  VHF  radio  pictures  of  lightning," 

J.  Geophys.  Res.,  76,  1478-1489,  1971. 


5.  M.J.  Master,  M.A.  Uman,  Y.T.  Lin,  and 
R.B.  Standler,  "Calculation  of  lightning  return 
stroke  electric  and  magnetic  fields  above 
ground,"  J.  Geophys.  Res.,  86,  12,127-12,132, 
1981. 

6.  Y.T.  Lin,  M.A.  Uman,  and  R.B.  Standler, 
"Lightning  return  stroke  models,"  J.  Geophys. 
Res.,  85,  1571-1583,  1980. 

7.  Y.T.  Lin,  M.A.  Uman,  J.A.  Tiller, 

R.D.  Brantley,  W.H.  Beasley,  E.P.  Krider,  and 
C.D.  Weidman,  "Characterization  of  lightning 
return  stroke  electric  and  magnetic  fields  from 
simultaneous  two-station  measurements,"  J. 
Geophys.  Res.,  84,  6307-6314,  1979. 

8.  M.A.  Uman,  C.E.  Swanberg,  J.A.  Tiller, 
Y.T.  Lin,  and  E.P.  Krider,  "Effects  of  200  km 
propagation  on  Florida  lightning  return  stroke 
electric  fields,"  Radio  Science,  11,  985-990, 
1976. 

9.  C.D.  Weidman,  and  E.P.  Krider,  "Sub¬ 
microsecond  risetimes  in  lightning  return-stroke 
fields,"  Geophys.  Res.  Lett.,  7,  No.  11,  955- 
958,  1980. 

10.  C.D.  Weidman,  and  E.P.  Krider,  "Sub¬ 
microsecond  risetimes  in  lightning  radiation 
fields,"  in  Lightning  Technology,  NASA  Confer¬ 
ence  Publication  2128,  FAA-RD-80-30,  1980. 


N 


Fig.  1  -  The  experimental  setup  including  the 
position  of  the  WC-130  at  17:09:40  EDT  on  July 
19,  1981  and  a  drawing  of  the  two  lightning 
channels  to  ground  deduced  from  VHF  time-of- 
arrival  measurements 
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Fiq.  1  -  Charqe  distribution  on 

charged  leader 

a.  Stepped  leader 
0.  After  returen  stroke 

c .  Cloud  di scharqe  advancinq  upwards 

d.  Cloud  discharge  adwandmc  downwards 


Fiq.  2  -  Charge  distribution  on 
uncharged  leader 

a.  Stepped  leader 

b.  After  return  stroke 

c.  Cloud  di  scharqe  beqirininq  stage 

d.  Cloud  discharge  end  stage 


Fig.  Z  -  Corona  discharge  on  Grbiter 


Fig.  4  -  Flash  over  on  Arbiter 
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Calculations  of  Lightning  Return  Stroke  Electric 
and  Magnetic  Fields  Above  Ground 
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The  first  detailed  calculations  of  lightning  return  stroke  electric  and  magnetic  fields  above  ground 
are  presented.  Waveforms  are  given  for  altitudes  from  0  to  10  km  and  ranges  trom  20  m  to  10  km 
These  waveforms  are  computed  using  the  model  of  Lin  el  al  ( 19801  and  a  modification  of  that  model  in 
which  the  initial  current  peak  decays  with  height  above  ground  both  the  original  and  the  modified 
models  result  in  accurate  prediction  of  measured  ground-based  fields  Rcium  siroke  held  measure¬ 
ments  above  ground  close  to  the  stroke,  with  which  the  calculations  could  be  compared,  have  nol  yet 
been  made.  Salient  aspects  of  the  calculated  fields  are  discussed,  including  their  use  in  calibrating 
airborne  field  measurements  from  simultaneous  ground  and  airborne  data 


Introduction 

Lin  el  al.  1 1 980J  have  recently  introduced  a  lightning 
return  stroke  model  with  which  return  stroke  electric  and 
magnetic  fields  measured  at  ground  level  I Lm  et  al..  1979) 
can  be  reproduced.  Here  we  use  that  model  and  a  modifica¬ 
tion  of  it  to  compute  electric  and  magnetic  fields  at  altitudes 
up  to  10  km  and  at  ranges  from  20  m  to  10  krn.  These 
calculations  provide  the  first  detailed  estimates  of  the  return 
stroke  fields  that  exist  above  ground  and  that  are  encoun¬ 
tered  by  aircraft  in  flight.  The  most  recent  generation  of 
aircraft  may  be  particularly  susceptible  to  lightning  electric 
and  magnetic  fields  because  these  aircraft  are  controlled 
with  low-voltage  digital  electronics  and  are  in  part  construct¬ 
ed  of  advanced  composite  materials  which  provide  a  reduced 
level  of  electromagnetic  shielding  (Corbin.  1979|.  Hence,  in 
the  context  of  aircraft  safety,  calculations  of  the  magnitude 
and  waveshape  of  airborne  electric  and  magnetic  fields  are  of 
considerable  practical  interest.  Furthermore,  from  a  scien¬ 
tific  point  of  view,  since  airborne  electric  and  magnetic  fields 
are  presently  being  measured  [Pius  ct  al..  1979.  Puts  and 
Thomas.  1981.  Baum.  |980|.  a  comparison  of  the  calcula 
Uons  given  in  this  paper  with  appropriate  experimental  data, 
when  they  are  available,  will  constitute  a  test  of  the  return 
stroke  model. 


Theory 

The  lightning  return  stroke  current  is  assumed  to  flow  in  a 
thin,  straight,  vertical  channel  of  height  H  above  a  perfectly 
conducting  ground  plane,  as  shown  in  Figure  I .  The  electric 
and  magnetic  fields  at  altitude  l  and  range  r  from  a  vertical 
dipole  of  length  dz'  at  height  c'  and  arbitrary  current  Hz' .  D 
are 
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where  1 1  >  and  (2)  are  expressed  in  cv  lindricai  coordinates.  en 
is  the  permittivity  and  a.,  the  pel meability  of  vacuum,  and  all 
geomelnc.il  (actors  ate  illustrated  in  Figure  I  Equations  ill 
and  (2)  are  obtained  in  a  straightforward  manner  using  an 
approach  similar  to  that  of  Uman  cl  al.  1 1975)  In  1 1 1  the  first 
and  fourth  terms  are  generally  called  the  electrostatic  field, 
the  second  and  fifth  the  induction  or  intermediate  field,  and 
the  remaining  two  terms  the  radiation  field.  In  (2)  (he  first 
term  is  the  induction  field  and  the  second,  the  radiation  field. 
The  effects  of  the  perfectly  conducting  ground  plane  on  the 
electric  and  magnetic  fields  due  Vo  the  source  dipole  are 
included  by  replacing  the  ground  plane  by  an  image  current 
dipole  at  distance  beneath  the  plane,  as  shown  in  Figure 
1  ( Stratton .  1941).  The  electric  and  magnetic  fields  due  to  the 
image  dipole  may  be  obtained  by  substituting  R;  for  R  and 
for  2 '  in  ( I >  and  (2)  above. 

In  this  paper  we  examine  only  thi  fields  of  a  typical 
subsequent  return  stroke  because  it  is  subsequent  strokes 
with  which  Lin  rt  al.  1 1980)  have  tested  their  model.  Subse¬ 
quent  strokes  are  more  easily  modeled  than  first  strokes 
because  m  contrast  to  firsts,  subsequenls  have  few  if  any 
branches,  have  relatively  constant  return  stroke  velocities, 
and  are  probably  initialed  at  ground  level  rather  than  by 
upward-going  leaders  ( Schonland ,  1956], 

The  model  of  Lin  et  al.  1 1980)  postulates  that  the  return 
stroke  current  is  composed  of  three  components:  (Da  short- 
duration  upward-propagating  pulse  of  current  of  constant 
magnitude  and  waveshape  associated  with  the  electrical 
breakdown  at  the  return  stroke  wavefront  and  responsible 
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for  the  return  stroke  peak  current;  the  pulse  is  assumed  to 
propagate  at  a  constant  velocity;  (2)  a  uniform  current  which 
may  already  be  flowing  (leader  current),  an  assumption  we 
use  in  this  paper,  or  it  may  start  to  flow  soon  after  the 
commencement  of  the  return  stroke;  and  (3)  a  ‘corona’ 
current  caused  by  the  downward  movement  of  the  charge 
initially  stored  in  the  corona  sheath  around  the  leader 
channel  and  discharged  by  the  passage  of  the  return  stroke 
wavefront.  These  three  current  components  are  illustrated  in 
Figure  2. 

Two  observations  form  the  basis  for  a  modification  of  the 
model  of  Lin  et  al.  1 1980]; 

1.  Al  the  lime  that  the  research  of  Lin  et  al.  11980]  was 
done,  subsequent  strokes  were  thought  to  have  both  lumi¬ 
nosities  (hence,  by  implication,  currents)  and  return  stroke 
velocities  that  were  invariant  with  height  [Schontand,  1956]. 
However,  Jordan  and  Uman  [1980]  have  since  shown  that 
subsequent  stroke  initial  peak  luminosity  varies  markedly 
with  height,  decreasing  to  half-value  in  less  than  i  km  above 
ground.  The  implication  of  this  result  is  that  the  breakdown 


pulse  current  (component  (1)  in  the  model)  must  also  de¬ 
crease  with  height. 

2.  When  the  breakdown  pulse  reaches  the  top  of  the 
channel,  the  model  of  Lin  el  al.  1 1980]  predicts  a  held  change 
of  opposite  polarity  to  that  of  the  initial  held,  the  waveshape 
of  the  field  change  being  a  ‘  mirror  image’  of  the  initial  field 
change.  A  detailed  discussion  of  the  mirror  image  effect  is 
given  by  Uman  ti  al.  ( 1975].  It  is  observed  occasionally  in 
the  fields  from  first  return  strokes  but  almost  never  in  the 
fields  from  subsequent  return  strokes  (Lm  ei  al..  1979],  If  the 
breakdown  pulse  current  is  allowed  to  decay  with  height  so 
that  it  has  a  negligible  value  when  it  reaches  the  end  of  the 
channel,  the  mirror  image  should  no  longer  be  manifest  in 
the  calculated  fields. 

In  view  of  observations  1  and  2  above,  we  propose  the 
following  modification  to  the  model  of  Lin  el  al.  1 1980]:  the 
breakdown  pulse  current  is  allowed  to  decrease  with  heigh; 
above  the  ground;  all  other  features  of  the  original  model 
remain  unchanged.  As  we  shall  see.  the  fields  at  ground  level 
produced  by  the  modified  model  are  essentially  the  same  as 
those  due  to  the  original  model  except  for  the  absence  of  the 
mirror  image.  However,  the  fields  in  the  air,  especial!)  at 
close  ranges,  differ  considerably. 

We  first  consider  the  calculation  of  the  electric  and 
magnetic  fields  of  a  typical  subsequent  stroke  using  the 
original  model  of  Lin  et  al.  |1980],  We  then  repeat  the 
calculation  for  the  modified  version  of  the  model.  The 
subsequent  stroke  used  in  this  study  is  that  for  which  the 
following  data  are  given  in  Figure  1 1  of  Lin  et  al.  1 1980):  both 
measured  and  calculated  fields  at  ranges  of  2  km  and  200  km 
at  ground  level  and  calculated  current  at  ground  level  This 
calculated  current  and  the  three  components  which  consti¬ 
tute  it  are  shown  in  Figure  3.  The  rise-to-peak  of  the 
breakdown  pulse  component  has  been  altered  from  that 
given  by  Lin  et  al.  1 1980]  so  as  to  be  consistent  with  the 
measurements  of  Weidman  and  Krider  11978],  The  salient 
parameters  of  the  current  used  in  the  field  calculations  for 
the  case  of  a  constant  breakdown  pulse  current  are  i  I 
breakdown  pulse  current:  increase  from  0  to  3  kA  in  1.0  as. 
followed  by  a  fast  transition  to  a  peak  value  of  14.9  kA  at  1.1 
ns.  half  value  at  3.8  ^.s,  and  zero  at  40  ms;  the  breakdown 


Fig.  2.  Current  distribution  for  the  model  of  Lin  et  at  1 1980]  in  which  the  breakdown  pulse  current  is  constant  with 
height.  The  constant  velocity  of  the  breakdown  pulse  current  is  v  Current  profiles  are  shown  at  four  different  times,  i, 
through  I,,  when  the  return  stroke  wavefront  and  the  breakdown  pulse  current  are  at  four  different  heights  ci  through  u. 
respectively. 
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Fig.  3.  Return  stroke  current  components  at  ground  calculated  from  measured  electric  and  magnetic  fields  for  a 

typical  subsequent  stroke 


pulse  current  propagates  upward  with  an  assumed  constant 
velocity  of  I  x  10®  m/s,  (2)  uniform  current:  3100  A,  and  (3) 
corona  current  tryected  per  meter  of  channel,  /oe 

with  l0  -  21  A/m,  \  =  1500m,  a  =  10-'  s'1,  and  b  = 
3  x  I06  s'1.  The  initial  charge  stored  on  the  leader  and 
lowered  by  the  return  stroke  is  0.3  C.  The  channel  length  is 
7.S  km.  In  the  modified  version  of  the  model  we  use  a 
breakdown  pulse  current  whose  amplitude  decreases  with 
height  as  e“"'\  with  k  ~  1500  m;  that  is,  the  breakdown 
pulse  current  decreases  with  height  in  exactly  the  same  way 
as  does  the  corona  current.  All  other  parameters  for  the 
modified  model  are  the  same  as  those  for  the  original  model. 

Since  both  first  and  subsequent  strokes  probably  have 
initial  currents  which  decrease  with  height  [Schonland.  1956; 
Jordan  and  ilman,  1980),  and  since  the  measured  wave¬ 
shapes  of  first  and  subsequent  stroke  fields  at  ground  level 
are  qualitatively  similar  f Lin  et  al..  1979],  one  would  also 
expect  the  airborne  subsequent  stroke  fields  calculated  using 
the  modified  model  to  be  qualitatively  similar  to  airborne 
first  stroke  fields. 

Results 

Calculated  vertical  and  horizontal  electric  fields  are  shown 
in  Figures  4 a  and  Ah.  respectively,  and  calculated  magnetic 
fields  in  Figure  4c.  Solid  lines  represent  the  original  model  of 
Lin  et  al.  1 1980)  and  dashed  lines  the  modified  version  of  the 
model  in  which  the  breakdown  pulse  current  decreases  with 
height.  All  zero  times  on  the  figures  indicate  the  time  at 
which  the  return  stroke  current  originates  at  ground  level. 
The  waveforms  at  the  field  points  begin  after  the  appropriate 
propagation  time  delay.  The  intersections  of  the  slanted  solid 
lines  with  the  horizontal  dotted  lines  at  various  heights 
indicate  the  times  at  which  the  return  stroke  wavefront 
passes  those  heights.  A  number  of  features  of  the  calculated 
waveforms  are  worthy  of  note: 

I .  With  the  exception  of  the  absence  of  the  abrupt  field 
changes  associated  with  the  end  of  the  channel,  the  fields  on 
the  ground  at  all  distances  and  the  fields  on  the  ground  and  in 
the  air  beyond  about  10  km  are  not  much  influenced  by  the 
breakdown  cunent  pulse  decrease  with  height.  This  is  the 
case  because  the  initial  parts  of  these  field  waveforms  are 
radiated  by  the  breakdown  current  pulse  while  it  is  very 
dose  to  the  ground  and  near  its  maximum  amplitude,  while 
later  portions  of  the  field  waveform  are  primarily  due  to  the 


uniform  and  corona  currents  which  are  the  same  in  the 
original  and  modified  models.  Hence  it  follows  that  ground- 
based  or  distant  airborne  measurements  cannot  be  used  to 
test  the  validity  of  the  model  modification  introduced  here. 

2.  The  abrupt  field  changes  associated  with  the  unatten- 
ualcd  breakdown  pulse  current  reaching  the  idealized  ends 
of  the  real  and  image  channels  \Uman  et  al..  1975]  do  not 
occur  when  that  pulse  current  is  attenuated  with  height  so 
that  it  does  not  have  an  appreciable  magnitude  when  it 
reaches  the  channel  end.  As  noted  earlier,  the  fact  that  these 
abrupt  changes  do  not  often  occur  in  the  experimental  data 
( Lin  et  al..  1979]  was  one  of  the  reasons  for  the  proposed 
modification  to  the  original  model  of  Lin  et  al  [1980). 

3.  At  ranges  less  than  about  200  m  the  horizontal  electric 
and  the  magnetic  field  components  above  ground  attain 
initial  peak  values  at  about  the  time  at  which  the  return 
stroke  breakdown  pulse  current  is  at  the  same  altitude  as  the 
field  point.  The  vertical  electric  field  component  undergoes  a 
sharp  decrease  at  this  point.  The  maximum  electnc  and 
magnetic  fields  are  due  to  the  charge  and  current,  respective¬ 
ly,  associated  with  the  breakdown  pulse  component  at  about 
the  lime  of  its  closest  approach  to  the  field  point.  The  peak 
electric  field  is  essentially  electrostatic,  the  peak  magnetic 
field  essentially  induction.  The  initial  peak  field  present  in 
measurements  made  beyond  a  few  kilometers  and  associated 
with  the  radiation  field  component  of  the  breakdown  pulse 
current  is  present  in  the  close  electric  and  magnetic  fields  but 
is  small  compared  with  the  electrostatic  and  induction  fields, 
respectively.  The  effect  of  the  decrease  of  the  breakdown 
pulse  current  with  height  is  primanly  to  decrease  the  magni¬ 
tude  of  the  initial  electrostatic  and  induction  peaks. 

4.  At  ranges  less  than  about  200  m  the  vertical  electric 
field  above  ground  is  bipolar  for  the  unattenuated  pulse 
current  due  to  the  passage  from  below  to  above  of  the  charge 
associated  with  the  breakdown  pulse  current.  As  one  moves 
farther  away  from  the  channel,  is  near  the  ground,  or 
considers  a  pulse  which  decays  with  height,  this  bipolar 
effect  is  reduced.  On  the  ground  near  the  channel  the  electric 
field  is  always  unipolar  and  of  opposite  polarity  as  compared 
to  the  initial  bipolar  field,  since  the  charge  motion  is  always 
at  a  height  above  that  of  the  field  point. 

5.  At  ranges  less  than  about  1  km  the  peak  value  of  the 
horizontal  electric  field  above  g round  is  larger  than  the 
associated  vertical  electric  field.  The  horizontal  and  vertical 
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Fig  4 a  Calculated  vertical  electric  fields  for  a  typical  subsequent  return  stroke  The  solid  lines  represent  the 
original  model  in  which  the  breakdown  pulse  current  is  constant  with  height,  the  dashed  lines  the  modified  modci  in 
which  it  decreases  with  height. 


fields  above  ground  are  roughly  equal  in  peak  magnitude  in 
the  3  km  range,  and  the  vertical  field  is  larger  beyond  about 
10  km. 

6.  At  ranges  greater  than  10  km  the  magnetic  field  and 
the  vertical  electric  field  are  relatively  weak  functions  of 
attitude,  whereas  the  horizontal  electric  field  increases 
roughly  linearly  with  altitude.  The  magnetic  field  and  the 
vertical  electric  field  are  height  independent  as  long  as  the 
difference  between  the  propagation  paths  from  the  source 
dipole  and  its  image  is  much  less  than  the  wavelength  of  the 
highest  significant  frequency  component  of  the  electromag¬ 
netic  radiation  from  the  source,  a  condition  which  is  met  to  a 
reasonable  approximation  at  a  range  of  about  10  km  at 
altitudes  below  about  3  km  for  the  current  waveshapes  used 
in  the  model.  Hence  measurements  of  distant  magnetic  and 
vertical  electric  fields  made  simultaneously  on  the  ground 
and  in  the  air  provide  a  simple  means  of  calibrating  the 
airborne  measurements. 

7.  The  initial  nonzero  value  associated  with  the  wave¬ 
forms.  which  can  be  clearly  seen  in  Figures  4a,  46.  and  4c  at, 
for  example,  10  km,  is  due  to  the  induction  field  from  the 
uniform  current  component  (component  2  in  the  model) 
assumed  to  exist  at  the  time  at  which  the  breakdown  pulse 
current  is  initiated  at  ground  level.  We  associate  the  uniform 
current  with  the  dan  leader  which  preceeds  the  return 
stroke.  The  electrostatic  field  value  at  the  time  of  the 


initiation  of  the  return  stroke  at  ground  due  to  the  current 
flow  prior  to  that  time  is  unknown  and  hence  is  not  included 
in  the  calculated  fields.  In  the  work  by  Lin  cl  al.  11980]  the 
initial  field  value  was  plotted  as  zero,  since  the  actual  value 
could  not  in  general  he  determined  from  the  type  of  measure¬ 
ments  made  by  Lin  cl  al.  1 197V). 

Disc  USSION 

In  this  paper  we  have  presented  the  first  detailed  estimates 
of  airborne  electric  and  magnetic  fields  due  to  lightning  We 
have  used  the  original  model  of  Lin  cl  al.  |I980|  and  also  a 
modified  version  based  on  observations  of  Lin  el  ai.  [1979] 
and  of  Jordan  and  Uman  [1980].  The  new  version  of  the 
model  (1)  results  in  fields  which  do  not  exhibit  abrupt 
changes  associated  with  the  breakdown  pulse  current  reach¬ 
ing  the  top  of  the  channel  and  (2)  can  be  expected  to  produce 
an  initial  luminosity  which  decreases  with  height  above  tne 
ground.  Though  the  individual  currents  which  define  the 
modified  model  are  not  unique  (see  discussion  by  Lin  et  ai. 
[1980]),  it  is  likely  that  the  total  current,  which  results  in 
accurate  prediction  of  ground-based  fields  and  is  consistent 
with  (I)  and  (2)  above,  predicts  airborne  fields  which  are  at 
least  qualitatively  correct. 

We  have  modeled  the  return  stroke  channel  as  a  straight 
vertical  antenna.  An  actual  return  stroke  channel  is  charac- 


Fig  4 h  Calculated  horizontal  electric  fields  for  a  typical  subsequent  return  stroke  The  solid  lines  represent  the 
original  model  in  which  the  breakdown  pulse  current  is  constant  with  height;  the  dashed  lines  the  modified  model  ill 
which  it  decreases  with  height. 


tended  by  tortuosity  on  a  scale  from  less  than  1  m  to  over  I 
km  |e.g..  Evans  and  Walker,  1963;  Hill,  1968].  Hill  |I969J 
and  LeVine  and  Meneghini  ( 1978],  using  simple  models, 
have  investigated  the  effects  of  channel  tortuosity  on  distant 
radiation  fields.  LeVine  and  Meneghini  find  that  the  wave¬ 
forms  computed  for  the  case  of  a  tortuous  channel  have  finer 
structure  than  those  for  a  straight  channel,  resulting  in  a 
frequency  spectrum  for  the  tortuous  channel  that  has  larger 
amplitude  at  frequencies  above  about  100  kHz.  Hill  shows 
that  horizontal  channel  sections  radiate  significantly  for 
frequencies  above  20-30  kHz  but  does  not  compare  the 
radiation  from  horizontal  channel  sections  with  that  from 
vertical  sections.  The  effect  of  using  a  tortuous  channel  to 
model  the  lightning  return  stroke  fields  at  close  range  has 
been  investigated  by  Pearlman  [1979],  again  using  a  simple 
model.  His  results  indicate  that  channel  tortuosity  has  little 
i  ifect  on  the  close  fields.  Since  the  peaks  in  the  close  electric 
and  magnetic  fields  are  due  to  the  charge  and  current, 
respectively,  associated  with  the  breakdown  pulse  current 
(as  discussed  in  (3)  of  the  previous  section),  we  suggest  on 
physical  grounds  that  the  general  shapes  of  the  close  fields 
should  not  be  greatly  different  from  those  shown.  However, 
the  peak  fields  at  close  range  should  occur  at  the  time  the 
breakdown  pulse  reaches  the  point  of  closest  approach  to  the 
held  point,  and  thus  the  distance  of  closest  approach  re¬ 
places  the  range  in  Figures  4a,  4b,  and  4c. 


Available  data  on  airborne  field  measurement--  are  limned 
to  the  observations  of  Pills  and  Thomas  |  1980)  and  of  Baum 
1 1980],  Pitts  and  Thomas  do  not  appear  to  have  any  data  on 
return  stroke  fields.  Baum  presents  airborne  measurements 
made  on  first  and  subsequent  return  stroke  electric  and 
magnetic  fields.  He  gives  one  typical  first  and  one  typical 
subsequent  return  stroke  electric  field  waveform  He  does 
not,  however,  make  an  independent  measurement  of  the 
distance  to  the  lightning  flashes  he  records  Rather,  he  uses 
the  values  of  the  observed  airborne  initial  peak  fields  and  the 
average  observed  values  on  the  ground  as  a  function  of 
distance  obtained  by  Lin  ei  at.  |  1979)  to  estimate  the  range. 
We  have  shown  in  this  paper  that  the  peak  radiation  fields 
for  distances  beyond  about  10  km  are  about  the  same  in  the 
air  and  on  the  ground.  However,  the  comparison  with 
average  values  of  the  fields  on  the  ground  as  a  function  of 
distance  can  lead  to  range  errors  of  two  or  three,  since 
individual  field  values  may  differ  from  the  average  by  this 
factor  [Lin  el  a!.,  1979].  If  we  do  use  Baum's  ranging 
technique,  the  subsequent  stroke  waveform  he  gives  is  at  a 
range  of  about  20  km.  The  aircraft  was  at  an  altitude  between 
3  km  and  5  km.  The  measured  airborne  electric  field  is  very 
similar  to  typical  measured  fields  on  the  ground  at  that 
range,  as  expected  from  the  theory  presented  in  this  paper. 
A  test  of  the  validity  of  the  predicted  fields  and  hence  the 
model  awaits  measurement  at  close  range  of  simultaneous 


12.132 


Master  et  al.:  Lightning  Return  Stroke  E  and  M  Fields 


Fig  4c.  Calculated  magnetic  fields  for  a  typical  subsequent  return  stroke  The  solid  lines  represent  the  original 
model  in  which  the  breakdown  pulse  current  is  constant  with  height;  the  dashed  lines  the  modified  model  in  which  it 
decreases  with  height. 


ground-based  and  airborne  return  stroke  electric  and  mag¬ 
netic  fields. 
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A  Comparison  of  Lightning  Electromagnetic  Fields  with 
the  Nuclear  Electromagnetic  Pulse  in  the  Frequency 

Range  104-10  Hz 
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Abstract — The  electromagnetic  fields  produced  h>  both  direct 
tiKhlnini*  strikes  and  nearbv  light  nine  are  compared  with  the  nuclear 
electromagnetic  pulse  iNKMPl  from  an  exnatmospheric  burst.  Model 
calculations  indicate  that,  in  the  frequence  range  from  III4  to  near  10 
Hz.  the  Fourier  amplitude  spectra  of  the  return  stroke  magnetic  fields 
near  ground  I  m  from  an  average  lightning  strike  will  exceed  that  of 
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the  Nf. Ml*.  Nearbv  first  rduin  strokes  at  j  range  of  about  Mi  in,  i! 
the\  are  severe,  product  electric -field  spectra  mar  groiimi  >  ■.  : 
exceed  that  of  the  NFMP  below  about  III'’  H/.  while  t h «  >pet;:. 
average  nearbv  first  return  strokes  exceed  that  of  tht  \fMf“  »*«  u.w 
about  3  *  10'  Hi.  Implications  of  these  results  | ■  »r  airiraft  m  P •  ? 

are  discussed. 

A>v  Words  —  l  ightning.  FMP.  electromagnetic  fields,  ampluiidi 
spectra,  aircraft. 


I  INTRODUCTION 

SERIES  ot  recent  articles  describe  the  threat  t,  the  l  r , . i. . 
States  command,  control,  and  communication'  it  1  i  : ,o • 
work  from  a  nuclear  electromagnetic  pulse  (I  MP  ,  \l  MI’i 
generated  by  an  cxoatinosphcrtc  nuclear  explosion.  Broad 
[l|-[3|  has  critically  examined  the  various  technical  an.:  > 
lineal  problems  surrounding  lerner  [4  J .  [5j  \i,  do 

cussed  the  damaging  crtccts  of  NfcMI’  on  the  C3  network  and 
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mulea:  power  plants:  and  Raloff  |6|,  (7|  has  considered 
.  : lucat  posed  by  NEMP  and  the  implementation  of  dcien- 
- : \ e  outeeies.  It  has  often  been  stated  |4|,  |6).  [S]  that  the 
i  lUMsat  NL  Ml’  are  comparable  to.  or  greatly  exceed,  those  of 
trie  most  severe  lightning.  For  example,  Holden  (8)  states  that 
"the  FMP  is  a  microsecond  burst  of  electromagnetic  energy,  a 
hundred  times  more  powerful  than  a  lightning  bolt."  As  far  as 
we  are  aware,  the  claims  that  NEMP  effects  almost  always  ex¬ 
ceed  those  of  lightning  arc  not  quantitatively  justified  in  the 
literature  On  the  other  hand,  a  recent  letter  to  the  Editor  of 
the  IEEE  Sm'THUM  |h]  cites  references  to  the  lightning  liter¬ 
ature  to  support  the  view  that  lightning  effects  can  be  equiva¬ 
lent  to  or  exceed  those  of  the  NEMP. 

The  effects  of  NEMP  lroni  an  exoatmospheric  burst  will  be 
toll  ovei  a  large  geographical  area,  whereas  the  effects  of  a 
single  lightning  discharge  are  local.  Nevertheless,  the  Ire- 
queries  ot  direct  and  nearby  strikes  to  sensitive  earthbound 
structures  like  nuclear  power  plants,  to  electric  transmission 
and  distribution  systems,  and  to  aircratt  in  flight  is  sufficiently 
high  ts>  warrant  a  careful  assessment  of  the  lightning  hazard 
Here  we  present  frequency-domain  comparisons  of  the  elec- 
i:i,  and  magnetic  fields  near  ground  due  to  model  lightning 
leiurn  strokes  with  those  of  the  NEMP  from  an  exoatmos- 
t'horic  burst  The  applicability  of  these  results  for  altitudes  at 
which,  lot  example,  aircraft  operate,  is  presently  a  matter  of 
'■peculation  due  to  the  paucity  of  airborne  measurements,  as 
w-  will  discuss  We  will  show  (bat,  in  the  frequency  range  from 
iO4  to  near  l(f  Hz,  the  calculated  Fourier  amplitude  spectra 

the  icttur:  stroke  magnetic  fields  near  ground  1  m  front  an 
.neruce  lightning  strike  wall  exceed  that  of  the  NEMP.  and  that 
c.e.-tiic  Held  spectra  neat  gionnd  ot  severe  nearby  first  return 
o  .  m  s  at  50  m  exceed  that  ot  the  NEMP  below  about  10**  Hz 
::k;  spectra  of  average  nearby  first  return  strokes  are  greater 
••clow  about  1  \  10*  Hz.  To  the  extent  that  fields  in  the  fre¬ 
quency  ramies  m  which  lightning  spectra  exceed  that  of  NEMP 
-cpiesent  a  hazard  b\ .  for  example,  exciting  resonances  m  a 
•»ru.  Hire  v.  Inch  couple  damaging  voltages  and  currents  to  dec- 
tioiiu  s  m  the  mtenor  of  that  structure  1 10| .  1 1  I ) .  lightning  ef- 
lects  can  apparently  be  as  severe  as  those  due  to  the  NEMP. 

II.  LIGHTNING 

Recently,  n  has  been  reported  that  the  electric  and  magnetic 
: te US  produced  by  all  important  lightning  discharge  processes 
contain  significant  variations  on  a  submicrosccond  time  scale 
[I2|-j  |x|  The  existence  of  these  field  components,  in  turn, 
implies  that  the  currents  which  produce  them  contain  large 
suhrmcroseconJ  variations  (Id),  (15).  ( !6J,  (|b ).  The  few 
wired  wide  band  measurements  of  lightning  currents  during 
sirixes  to  airplanes  in  (light  show  submicrosccond  rise  times 
tor  dirrcnt  pulses  in  the  100- A  range  [20[ ,  |21 ).  These  pulses 
are  probably  associated  with  small  cloud  discharge  processes. 

I  man  et  al  [22)  have  calculated  the  distant  electric  and 
magnetic  radiation  fields  produced  at  ground  level  by  a  fixed 
current  waveshape  propagating  up  a  straight  vertical  channel 


(D,  r)  = 

•  Hr  -  D/c)at, 

t^D/c 

(1) 

2nD 

V 

t  i*  D/c 

,( D .  t)  --- 

- •  i(r  -  D/c)a0, 

2nd) 

(2) 

where  r(Z),  t  ?-■  0,  is  (he  current  at  lime  /.  v  is  the  velocity  w  ill, 
which  the  current  pulse  propagates  up  the  channel.  l>  i\  the 
horizontal  distance  from  the  channel  to  the  point  at  wincii  tin 
field  is  measured,  c  is  i lie  speed  of  light.  p„  is  the  pe.-mcabii.ti 
ot  free  space,  r  is  the  vertical  coordinate,  and  o  is  the  azimuthal 
coordinate.  The  best  available  models  for  the  currcni  m  the  re¬ 
turn  stroke  phase  of  a  cloud- to-ground  discharge  1 2  > J  .  (24j 
(sec  [25],  [2bj  lor  a  general  review  of  lightning  dischaige 
phenomena)  have  model  currents  which,  in  the  time  domain, 
produce  electric  and  magnetic  fields  in  good  agreement  with 
wide  band  (dc  to  about  2  MHz)  time-domain  measuieinents 
made  at  ground  level.  For  these  models.  ( 1 1  and  (2)  piovide  a 
good  approximation  to  the  relation  between  the  initial  return 
stroke  radiation  field  and  the  initial  return  stroke  cunenr 
Weidman  and  lkndei  |)(>]  have  measured  the  maximum  rate- 
of-rise  of  the  initial  return  stroke  radiation  field  V:  rest 
strokes  and  find  a  mean  of  about  50  Von-ps  normalize,' 
an  inverse  range  relation  to  a  distance  of  I  (Hi  km.  This  u.o.e 
value,  when  substituted  m  (1).  with  an  assumed  return  sir., no 
velocity  of  10s  m  s.  leads  to  a  calculated  mean  maximum 
rate-of-rise  of  the  return  stroke  current  oi  ahoui  1'0  s\  ps.  a 
value  which  is  representative  of  the  current  iust  above  groin:  . 
The  maximum  values  of  maximum  ratc-of-rise  of  fielu  .mo  ct.t- 
rent  from  CHT  measured  iirst  strokes  are  about  2.5  times  :r.c 
mean  [  1 6 1 . 

Lightning  return  stroke  current  wjvctortm  nave  iu-ei.  *:t- 
rectly  measured  ilurinp  strikes  to  instrumented  imvc;-.  jo 
Swit/erland  [27j.  m  Italy  { 2  S ) .  and  m  South  An  tea  [2  ' 
I'lifotuinately.  cutrents  to  tall  towers  do  not  tnvessaiily  p;>- 
vide  a  good  estimate  oi  the  current  encountered  by  sr.ia’:. 
earthbound  structures  or  of  the  current  m  the  lightning  ck.ni 
nel  above  ground  necause  ot  the  effects  ot  the  relatively  Ion; 
upward-propagatmg  leader  which  is  initiated  by  the  tower  am: 
because  of  the  effects  of  the  tower  inductance,  capacttuitc- . 
and  relatively  large  ground  impedance  characteristic  o;  -.iu 
mountain,.-.. us  terrain  where  most  measurements  have  me: 
made.  The  upward-propaga i ini:  discharge,  lor  example,  c.  a:  1,1 
cause  a  slower  oveiall  cunent  use  time  at  the  towei  than  a 
comparable  strike  to  normal  ground  and  could  conceivably  al¬ 
ter  or  mask  the  fast  current  components.  This  effect  should  be 
more  pronounced  lot  lust  u-turn  strokes  than  to:  subsequent 
strokes,  since  the  latte:  ate  thought  not  to  have  long  upward- 
propagating  leaders.  Uergei  ft  al.  [27]  found  that  me  median 
peak  current  foi  first  return  strokes  which  lowered  negative 
charge  to  a  tower  m  Switzerland  was  50  kA  and  that  the 
median  maximum  raie-of-n.se  of  the  .-ttneni  was  12  kA  us. 
The  corresponding  values  at  the  5. percent  level  were  SO  kA 
and  52  kA/jus.  For  negative  stiokes  subsequent  to  the  first 
multiple-stroke  Hashes,  the  median  peak  current  was  12  kA. 
and  the  median  maximum  ratc-ot-nse  was  40  kA  us.  Subse¬ 
quent  strokes  at  t he  5-pcrccnt  level  bad  a  peak  current  of  50 
kA  and  a  maximum  ratc-of-rise  of  120  kA/ps.  It  is  interesting 
to  note  that  the  subsequent  stroke  currents  reported  in  [27] 
have  shorter  overall  rise  times  and  larger  maximum  raies-of- 
rtse  than  first  strokes.  This  result  should  be  contrasted  with 
the  electric  radiation  field  measurements  made  by  Weidman 
and  Krider  [13|,  (16),  who  report  no  significant  differences 
in  the  maximum  rates-of-rise  of  first  and  subsequent  stroke 
fields.  The  tower  on  which  Berger  <■;  al.  [27]  nude  their 
measurements  was  on  top  of  Mi.  San  Sjlvatoie  ir.  Switzerland 
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TABLE I  A 

Current  Parameters  for  an  average  first  return 
STROKE  OBTAINED  FROM  REMOTE  MELD  MEASUREMENTS 
IN  ACCORDANCE  WITH  THE  PROCEDURE  OUTLINED  IN 
1231.124) 


Curr«nt  at  ground 


Time  (  ws ) 


Current  tkA) 


0.0 

0.0 

100.  0 

5.0 

105.0 

20.0 

105.1 

35.0 

107.  o 

18.o 

1»2.5 

25.0 

uo.o 

27.0 

llil.l 

18.0 

ihij.n 

12  .*• 

2OO.0 

5.0 

30o.o 

0.0 

gleen  hy  the  to  I  Low  in/,  paraattcrii 
(1)  Breakdown  puls*  current  with  velocity  1 
t  (  us ) 

0.  i» 

5.0 

S.l 

I'l.ft 

•4ll.lt 


The  pulse  decay*  exponentially  wtin  height  above  the  ground 
with  a  decav  constant  •  2.')  x  10*  m 

(2)  Irutorm  cur  lent  «u  •  5  kA 

time  duration  -  ‘.3  m;  turn—*n  time  ■  U  is 

(1)  t-oron*  r.rrent  prr  unit  lenpth  is  Jf  -  _  r~z 

where , 


tabu:  ii  a 

C  urrent Parameters  tor  an  average  subs!  yui  nt 
Return  stroke  obtained  from  Rfmoti  i  ield 

MEASUREMENTS  |2.t|.  [24) 


Current  at  ground 

Time  ( *.&  i 

100.  o 

J  i ;  1 .  > 

l  '3.8 
1  »b .  5 

12.'.  n 

2- >0.0 
3  iV:.  • 

Riven  bv  the  following  parameters 

(1)  Breakdown  poise  current  with  veloci: 

t ;  -s) 


; ;  .  r- 

The  pulse  decavs  exponentially  with  height 
with  a  decay  constant,  .  5  x  1  '  e, 

'  '  Pnttorm  current  »  3  k\ 

ttts-  duratio,,  -  ■. 

i 3)  Corona  cut  rent  per  uni l  lenct' 
where  , 


Current  . <A 

3.1 


t  urn -on  tio 


»  S->.  'i  V  m 
4  *  ?,’»  <  in 4  m 


U  '  l  A  I  '  N*- 

?  -  1  x  l  '"  s  ' 


TABLE IH 

Current  Parameters  for  an  average itrst  return 

STROKE  OBTAINED  FROM  TOWER  MEASUREMENTS  |27| 


T  l  me  v  ws  ) 

Current 

0.0 

<».•’ 

3.0 

3.2 

5.0 

7.0 

8.5 

14.0 

7.8 

21.o 

8.3 

?« .  0 

«.* 

32.0 

9.0 

34.0 

10.0 

35.0 

40.0 

2fl.  n 

SO." 

21.*) 

25U.O 

•J.o 

On  the  other  hand,  the  tower  used  in  the  South  African  study 
(29]  was  situated  on  a  comparatively  flat  terrain.  In  (29],  a 
maximum  current  rate-of-rise  of  180  kA/gis  was  reported  for  a 
subsequent  stroke,  although  the  total  sample  size  was  only  1 1 
flashes. 

Berger  ei  al.  [27j  have  computed  “mean  lightning  current 
waveshapes"  for  first  and  subsequent  negative  strokes.  These 
were  determined  by  first  normalizing  all  the  measured  wave- 


TAHLI  II  B 

Current  Parameters i ok  an  avi rage  si hsi  yt 

RETURN  STROKE  OBTAINED  FROM  TOWI  i 
MV  ASU REM ENTS  [2?  | 


I.urr 


1- 

1  fl.n 
1 


forms  to  a  peak  value  of  unity  and  then  averaging  the  meas¬ 
ured  values  at  each  tune.  If  this  "mean  waveshape"  is  scaled 
up  in  current  to  produce  a  large-amplitude  waveform,  as  we 
shall  do  to  model  a  severe  lightning,  the  rate-of-nse  necessarily 
scales  also.  On  the  other  hand,  the  reported  tower  measure¬ 
ments  do  not  show  a  very  strong  correlation  between  the  peak 
current  and  the  rate-of-rise  of  current  (27( ,  (30] . 

There  is  general  agreement  that  the  mean  peak  current  dur¬ 
ing  strikes  to  normal  objects  on  the  ground  is  in  the  20-40-kA 
range,  and  that  peak  currents  of  175-.  25  kA  are  present  in 
about  1  percent  of  the  strikes  (31] . 
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1  o  calculate  return  stroke  fields  tor  comparison  with  the 
M  Ml’,  we  will  use  first  and  subsequent  return  stroke  current 
wavelorms  that  are  interred  from  both  the  remote  electro- 
magnetic-field  measurements  and  the  tower  measurements. 
Specification  of  currents  for  average  first  and  subsequent 
strokes  is  given  in  Tables  I  and  II.  The  peak  current  value  for 
an  average  first  stroke  is  chosen  to  be  35  kA,  and  for  an  aver¬ 
age  subsequent  stroke,  18  kA.  The  currents  derived  from  the 
electromagnetic-field  measurements  |23],  [24],  differ  from 
the  directly  measured  currents  on  Ml.  San  Salvatore  [27]  pri¬ 
marily  in  the  relatively  slow  rate-of-risc  of  first  stroke  current 
in  the  tower  measurements  compared  to  that  derived  from  the 
fields.  Severe  lightning  currents  are  obtained  by  multiplying  by 
a  factor  of  jive  both  the  typical  currents  derived  from  electro¬ 
magnetic  fields  |23|,  |24|  and  the  "mean  lightning  current 
waveshapes"  trom  tower  measurements  [27|  given  in  Tables  I 
and  It  Although  the  peak  value  of  the  currents  determined 
this  way  are  representative  of  measured  severe  liglunuig  [31], 
the  rate-of-rise  we  use  lor  severe  lightning,  live  times  the  aver¬ 
age  value,  may  be  excessive  it  the  rate-of-change  of  current 
does  not  scale  with  the  current.  As  noted  previously  .  no  strong 
correlation  has  been  found  between  peak  current  and  rate-of- 
change  of  current  in  the  tower  measurements  J 2” ]  ,  |30j.  No 
uata  have  been  published  on  this  correlation  for  the  currents 
derived  from  the  fields,  and.  as  noted  previously,  the  largest 
value  ot  the  maximum  ruie-of-rise  of  the  electric  radiation  field 
tor  4"  jirst  strokes  was  only  about  2.5  times  the  mean  |  lo| . 

III.  NFMP 

The  characteristics  of  the  NFMP  are  a  function  ot  whclhet 
'.lie  iiucie.ii  event  is  in  ot  out  ot  the  atmosphere  and  the  dis¬ 
tance  trom  it.  A  thorough  survey  ot  the  mechanisms  by  which 
tne  NfcMI’  is  geneiated,  the  details  ol  the  coupling  of  the 
NT-' Ml’  to  a  variety  ol  sy  stems,  and  the  response  of  those  sys¬ 
tems  is  found  in  two  volumes  of  collected  papers  [10],  |l  ij 
Reasonable  approximations  to  the  NFMP  waveform  at  the  sui- 
faee  of  the  earth  ot  at  aircraft  altitude  due  to  an  exoatmos- 
phenc  burst  have  been  given  in  I.ee  ( I  I  | .  For  the  present  study , 
we  choose  the  exoatmospherie  burst  NFMP  waveform  trom 
Lee  [  1  1 1  which  appears  to  be  the  choice  ot  most  NT  MP  re¬ 
searchers 

],  t>0  (3) 

//(/)  --  //0[t  ‘u\ ,  i  0  t4) 

with  /: n  =  5.2  X  I04  V/m,  ll0  =  1.4  X  102  A /V  -  4.0  X 
10°  s  1 .  and  d  -  5.0  X  10h  s"  '. 

IV.  COMPARISONS 
A  Direct  Strike  and  MEM P 

To  compare  the  fields  from  lightning  direct  strikes  with 
NFMP  fields,  we  must  choose  an  example  object  to  be  struck. 
Let  us  consider  the  fields  at  the  surface  of  a  hypothetical 
cylindrical  metallic  aircraft  fuselage  of  radius  r.  Wc  choose  the 
atrcrjft  as  an  example  because  of  its  considerable  practical  im¬ 
portance.  If  a  lightning  return  stroke  attaches  directly  to  this 
aircraft  anJ  the  current  flows  uniformly  along  the  fuselage,  the 


magnetic-field  intensity  at  the  surface  will  be  about 
/ 

//  =  —  t  5 1 

2nr 

where  /  is  the  lightning  current  and  where  the  total  field  has 
been  approximated  as  magnetostatic.  Obviously,  the  field  will 
be  the  same  at  a  distance  r  from  the  axis  of  any  structure 
much  longer  than  r  which  uniformly  carries  the  current,  and 
hence. the  results  to  be  obtained  are  generally  applicable.  1: 
the  aircraft  is  struck  by  an  average  first  return  stroke  with  a 
peak  current  of  about  35  kA,  the  peak  magnetic  field  a', 
surface  with  an  assumed  radius  of  1  m  will  be  about  -  c.  ,■ 
103  A/m.  A  175-kA  severe  stroke  will  produce  a  peak  field  of 
about  2.x  X  104  A/in.  If  about  half  of  the  lightning  field  rises 
in  peak  in  about  0.1  s.  as  suggested  by  the  electromagnet',. - 
field  measurements  ofWcidman  and  Under  |!3],  [  1  (> ] .  tr.en 
the  maximum  late  of  change  of  ihe  magnetic  field  trom  at; 
average  first  stroke  will  be  about  2.8  X  1010  A.'nrsand  ::om 
a  severe  stroke  will  he  about  1.4  X  1 0 1  :  A.m-s  The  pus 
NFMP  field,  2.x  >  10:  A  m.  can  be  obtained  by  multiplyim 
(4)  by  i  factor  of  2.  to  take  account  ot  the  reflection  ot  \:.j 
NEMP  plane  wave  from  the  metallic  surface  of  the  attcra  ; 
The  maximum  field  late-ut-change  tor  NFMP  is  1.4  • 

'  m -s  and  exists  for  a  time  ol  the  older  ot  a  nanosecond.  'I 
return  stroke  peak  fields  trom  normal  lightning  exceec.  . 
of  the  NEMP  by  a  factor  of  about  20.  The  NFMP  maximum 
rate-of-change  exceeds  tiiat  of  normal  lightning  in  a  tactur 
about  5  and  is  about  equal  to  that  of  severe  lightning. 

We  now  examine  how  the  lime- domain  natamcteis  ce:;\e„ 
above  tor  lightning  and  NFMP  are  retlccted  in  tin  1  .  • 

am  nil  time  speed  a  tor  the  two  events.  A  •cam  w  e  use  ti.e  exaur.v 
ol  the  fields  on  the  surface  of  an  aircuti.  the  caicuiated 
nine  ticKis.  however,  being  essentially  the  same  jt  comparable 
distance  trom  ans  direct  strike.  For  the  comt’utufio::-  : 
volving  currents  derived  from  elect  tic  and  magnetic  ;k.: 
boll;  the  radiation  and  induction  field  terms  m  the  eer..m. 
field  equations  Have  beer,  included  ]22'].  ]24],  although  t5i. 
with  the  currents  given  m  Tables  IA  and  1 1 A .  provides  a  eonj 
approximation  to  the  end  result  In  calcui.it i. ns  invoivu,  •  cut- 
rents  trom  tower  measurements,  the  magnetic  lielu  is  can  u 
latcd  directly  from  (5). 

Figs.  I  and  2  show  the  Fourier  amplitude  spectra  o:  the 
time-domain  magnetic  fields  produced  tu  currents  given  ,r. 
Tables  1  and  II  for  both  average  and  severe  return  strokes  de¬ 
rived  from  both  tower  and  remote  field  measurements.  Tne 
current  waveforms  are  composed  of  straight-line  segments  be¬ 
tween  the  points  given  in  Tables  1  and  II  and  are  digifi/v,..  ; 
0.005-qs  intervals  for  the  calculation  ot  the  Fourier  amplitude 
spectra.  The  spectra  inferred  from  the  electromagnet!. -field 
measurements  are  larger  above  10'  H/  than  those  derived 
from  the  lower  data  for  first  strokes,  due  to  the  relatively 
slow  first-stroke  current  rate-of-rise  measured  on  towers,  hut 
the  two  spectra  are  similar  for  subsequent  strokes.  For  aver¬ 
age  return  strokes,  the  spectral  amplitudes  for  the  first  and 
subsequent  stroke  fields  determined  from  remote  electromag¬ 
netic  measurements,  and  the  subsequent-stroke  measured 
tower  current,  arc  equal  to  the  NEMP  at  a  frequency  near 
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Fig.  1  Magnetic-field  Fourier  amplitude  spectra  for  a  direct  strike  by 
an  average  and  a  severe  first  return  stroke  from  both  tower  and  re¬ 
mote  lield  measurements  and  for  NEMF.  The  dips  in  the  remote 
field  data  at  1  x  ID7  and  2  X  107  Hi  are  due  to  the  large  linear  cur 
rent  transition  taking  place  in  0  1  us  during  the  current  rise  to  peak. 
Since  real  waveshapes  do  not  have  such  linear  transitions,  the  dips 
in  the  spectra  arc  artificial. 

•40, 


>0'  't>4  •!>'  »0*  lO’  '0* 

HUOU€NCY  IHfl 

Fig.  2.  Magnetic-field  Fourier  amplitude  spectra  for  a  direct  strike 
by  an  average  and  a  severe  subsequent  return  stroke  from  both 
tow  er  and  remote  Field  measurements  and  for  NEMP.  The  dips  in  the 
remote  field  data  at  1  x  107  and  2  x  107  Hz  are  due  to  the  large 
linear  current  transition  taking  place  in  0.1  ms  during  the  current 
rise  to  peak.  Since  real  waveshapes  do  not  have  such  linear  transi¬ 
tions,  the  dips  in  the  spectra  are  artificial. 

107  Hz  and  exceed  the  NEMP  at  frequencies  below  that  value. 
For  above-average  return  strokes,  the  lightning  spectra  exceed 
the  NEMP  spectra  to  frequencies  above  107  Hz. 

The  preceding  comparisons  between  lightning  and  NEMP 
Fields  are,  strictly  speaking,  applicable  only  near  ground.  We 
use  the  example  of  an  aircraft  because  of  its  practical  impor¬ 
tance,  and,  for  that  reason,  a  discussion  of  the  Fields  above 
ground  is  in  order.  An  aircraft  in  flight  would  probably  not 
encounter  the  full  return  stroke  current  which  would  flow 
through  a  structure  on  the  ground  since  the  return  stroke  cur¬ 
rent  will  probably  decrease  with  height  [24],  In  fact,  Clifford 
and  Kasemir  [32]  argue  that  most  strikes  to  aircraft  are  not 
return  strokes  but  are  triggered  by  the  aircraft  and  are  some 
sort  of  in-cloud  discharge  with  a  rate-of-change  of  current  an 
order  of  magnitude  less  than  that  of  the  average  return  strokes 
we  are  considering.  Clifford  and  Kasemir  [32]  contend  that 
aircraft  are  only  occasionally  involved  with  return  strokes,  al¬ 
though  the  total  data  from  instrumented  aircraft,  on  which 
this  opinion  is  based,  are  meager.  In  any  event,  relatively  large 
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return  stroke  currents  at  ground  level  may  well  still  pro^.c 
currents  at  aircraft  operational  altitudes  to  cause  lielJs  equti- 
alent  to  the  NEMP  at  frequencies  below  about  107  Hz. 
important  is  the  observation  that  both  the  in-tloud  drsena;  cl- 
processes  which  precede  stepped  leaders  in  ground  flashes  aim 
certain  pulses  in  intracloud  lightning  discharges  procure 
Fourier  amplitude  spectra  measured  near  ground  lor  distant 
discharges  comparable  to  those  of  distant  return  strokes  ( 1 '  j . 
implying  that  there  are  in-cloud  events  which  produce  .lose 
fields  in  the  cloud  equivalent  to  close  return  stroke  fields  t  ea: 
ground.  These  in-cloud  processes  can  be  expected  i"  mteta.: 
with  aircraft.  An  accurate  assessment  of  the  probability 
craft  involvement  with  different  types  and  phases  ut  ;.cr, . 
awaits  further  research.  Finally,  it  is  worth  noting  tn.it  \ 
NEMP  wavefront  is  plane  while  the  hghtntnc  fient  j% 
and  that  lightning  channel  attachment  to  ar>  aircratt  nm  -lic- 
the  behavior  of  traveling  and  reflected  waves  on  the  j,:,:.'. 
structure  from  the  free  field  NEMP  case,  anu  hence  the-c 
be  additional  factors  in  the  comparison  which  we  !...ve  :n< 
considered. 

B.  .Xearbv  Lightning  and  XL'MP 

For  the  direct  lightning  strike,  we  have  compared  :1  to..: 

netre  field  at  an  aircraft  surface  with  the  NEMP.  Fo:  a 
flash,  we  will  compare  the  electric  fields.  The  fields  ate 
which  would  exist  in  the  absence  of  the  aircraft  Wt  ..,11  p. 
spectra  only  for  the  severe  first  return  stroke  at  c;  --:,.: 

In  Fig.  3,  we  show  the  NEMP  spectrum  calculated 
the  expression  given  in  (?)  along  with  three  electric  fici 
amplitude  spectra  for  severe  first  return  strokes  wt.s. 
strike  the  ground  50  m  from  the  observation  point.  The  ’luce 
lightning  amplitude  spectra  arc.  !)  the  average  first 
electric  radiation  field  spectrum  mcjsureu  by  U'ewrr.a;.  ' 
[17]  for  return  strokes  at  about  50  km  extrapoia'e-a  t 
using  an  inverse  distance  relationship  and  multiplier  ... 
tor  of  5  to  simulate  a  severe  stroke;  2)  the  clcetnc  ;acu:;.-:! 
field  spectrum  at  50  km  calculated  using  the  model  of  V. ,  ,;et 
et  al.  [24]  with  the  currents  in  Table  1A  multiplied  by  a  fac¬ 
tor  of  5  and  extrapolated  to  50  m  using  an  inverse  distance 
relationship;  3nd  3)  the  total  electric-field  spectrum  at  50  m 
calculated  using  the  model  [24]  with  the  currents  m  Tame 
LA  multiplied  by  a  factor  of  5. 

The  calculated  and  the  measured  radiation  ficid  spec;:*  at 
50  km  extrapolated  to  50  m  are  essentially  identical.  The  amp¬ 
litude  spectrum  computed  for  the  total  electric  field  of  a 
1 75-kA  lightning  at  50  m  is  equal  to  the  extrapolated  rad;*  ,.>n 
field  near  107  Hz  and  is  greater  for  lower  frequencies  u’causc 
Ihe  elecirastalic  and  induction  components  of  the  total  field 
add  to  the  radiation  component.  The  spectrum  of  the  total 
electric  field  exceeds  that  of  the  NEMP  below  about  10‘>  Hz. 
For  an  average  nearby  first  return  stroke,  the  total  electrical 
field  spectrum  exceeds  that  of  the  NEMP  below  about  3  X 
10s  Hz. 

The  nearby  discharge  has  been  assumed  to  be  at  a  distance 
of  50  m  because  this  is  about  the  range  at  which  an  earth- 
bound  structure  or  an  aircraft  in  flight  would  be  expected  to 
become  involved  in  a  typical  direct  strike.  The  closest  distance 
at  which  the  lightning  return  stroke  fails  to  attach  to  a  ground- 
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fig  3.  EleciruM’ield  Fourier  amplitude  spectra  lor  a  severe  first  re¬ 
turn  stroke  at  50  m  and  for  NEMP.  For  the  purpose  of  clarity,  dips 
in  the  computed  50-km  spectrum,  similar  to  those  shown  in  Figs. 
1  and  2.  have  been  omitted 


based  structure  tnav  be  derived  from  the  “striking  distance" 
concept  developed  by  Gold  [33 1 . 

Fig.  3  shows  the  vertical  electric  field  at  ground  level. 
Master  el  al.  [24)  have  computed  the  horizontal  and  vertical 
electric  fields  and  the  horizontal  magnetic  fields  above  ground 
level.  The  validity  of  these  calculations,  however,  depends  on 
the  validity  of  the  assumed  currents  above  ground  level.  Field 
measurements  above  ground  are  needed  before  the  adequacy 
of  the  calculated  nearby  return  stroke  field  environment  above 
ground  can  be  evaluated.  Further,  future  calculations  should 
include  the  effects  of  channel  tortuosity  which  arc  not  taken 
into  account  in  the  model  1 23 [ ,  [24)  used  in  this  paper.  Al¬ 
though  nearby  return  stroke  fields  may  well  decrease  with 
altitude  to  the  point  that  they  are  unimportant  compared  to 
;iie  NEMP.  in  cloud  processes,  as  noted  earlier,  can  produce 
distant  fields  near  ground  with  comparable  Fourier  amplitude 
spectra  to  return  strokes  [17],  implying  that  nearby  in-cloud 
processes  at  aircraft  altitudes  can  produce  a  similar  electro¬ 
magnetic  environment  to  that  of  nearby  return  strokes  at 
ground  level. 

V.  DISCUSSION  AND  CONCLUSIONS 

In  this  paper,  we  have  compared  calculations  of  the  electro¬ 
magnetic  environment  in  the  frequency  range  104-10'  Hz  pro¬ 
duced  near  ground  by  direct  and  nearby  lightning  return  strokes 
with  those  due  to  the  NEMP  produced  by  an  exoatmospheric 
burst.  In  that  frequency  range,  the  calculated  Fourier  amplitude 
spectra  of  the  magnetic  fields  1  m  from  a  direct  lightning  strike 
by  an  average  return  stroke  are  greater  than  that  of  the  NEMP. 
The  spectra  of  severe  nearby  first  return  strokes  at  about  50  m 
exceed  that  of  the  NEMP  below  about  10*  Hz  of  average 
nearby  first  strokes  below  about  3  X  10s  Hz.  These  results  fol¬ 
low  primarily  from  the  time-domain  return  stroke  currents 
given  in  Tables  IA  and  IIA,  which  were  determined  by  theory 
from  measured  time-domain  electric  and  magnetic  fields.  The 
frequency  spectra  calculated  from  fields  computed  using  these 
currents  are  in  excellent  agreement  with  measured  frequency 
spectra  up  to  107  Hz  for  lightning  at  a  distance  of  50  km,  as  is 
evident  from  Fig.  3,  providing  confidence  in  the  model  currents. 

Although  many  of  the  computations  in  this  paper  were  per¬ 
formed  for  an  idealized  aircraft  near  ground,  the  basic  results 


are  applicable  to  any  similar  size  ground-based  system  and,  to 
the  extent  discussed,  to  aircraft  at  flight  altitudes. 

REFERENCES 

[II  \5  J  Broad.  Nuclear  pulse  (li  Awakening  to  ihc  chaos  factor.' 
Set  .  v u I  212.  pp  1009-1012.  1981 

[2]  VV  J  Broad.  Nuclear  pulse  (II)  Ensuring  delivers  til  the  dooms¬ 
day  signal.  Sci  .  vol  212,  pp  1 1  16-1  1 20.  198  I 

(3)  W  J  Broad.  "Nuclear  pulse  (III)  Playing  a  wild  card.  '  5c/  .vol 
212.  pp  1248-1251 .  1981 

1 4  J  P  J  Lerner.  "Electromagnetic  pulses  Poiential  cripple?.  I  LEE 
Spectrum,  vol  IK.  pp  4  1—46  May  I  9K  I 

|5J  E  J  Lerner.  ‘  EMPs  and  nuclear  power.  /Lit  Span. nr  v.*i 
IK.  pp  4K— 49 .  June  1981 

[6j  J  Raloll  "EMP  A  sleeping  electrons  dragon.  V,  v,-. 

I  19.  pp  (no-  3t)2.  iwn  [ 

J  Rjlofl.  F.MP  Dctensive  strategic-  .Si.  v«>:  .  v.  pf 

314-  3 15  1 98  t 

J8|  C  Holden.  Energy,  security,  and  wai.  Si .  v ,  > ,  21  i  p  h.vv 
I9KI 

|9  j  P  B  Corn  and  J  C  Corbin.  Letter  t«'  the  editor.  It  It  >;n . 

irum.  vol  IK.  p  20  (At  19K1 

[10]  Special  Issue  ot  //  Li  t  runs  1  tromavn  t  omiust  voi  l  All  2(l 
Feb  |97k  (also  published  as  ILEL  Trans  Antenrius  Rropuvui  v... 
AP-2b.  Jan  19?x>. 

[Ilj  K  H  S  Lee.  Id  F.MP  interaction  Piuiviplc-  ic* rminue-  ai  •. 
reference  data.  LM/'  Interaction  2  i  \ f  V\  1  Ik  so  4;C  i>t  . 

1 980 

[12]  l  P  K rider.  (  P  Ske.Jma-..  and  R  C  \.Vi,e  lieu  .m:  , 

fields  produced  ov  lightning  stepped  leader-  ./  A. 

v  i*|  82  .  p['  9S  I  960  1 9"“ 

|  I  3)  C  D  Wcidnun  and  I  P  Node:.  The  one  strut,  ure 

return  stroke  waveforms  J  Geophw  Res  \o-  Kv  pp 
024'.  1 9  As 

[14]  C  [)  W  eidman  and  l  P  krtdot  Submit, i .‘second  rise  tores  :. 

produced  by  miracioud  lightning  discharge  processes  ./  (,c>- 
ph\s  Res  voi  s4.  pp  3IM  |9-g 

[  1 5 {  (  l)  U  ndiiun  and  F  /’  k  rider  Sunrti:d’< 'second  rise  .  n 

lightning  radiation  fields,  m  iiehimne  lecnnoioev  #i 
It'ih  Ssrnp  \AY3  tangles  Rrsetireh  i  ente' .  •  ••  >  N  AS  >. 

CP2I2K.  I  AA-R1)  KO  (0  Idso.  pp  :u  A 
j  1 6)  C  1)  Weidm.m  and  l  1’  kr.de;  Suomi*  t>.i  : . - u  ;r:  «•-  • 

iightmng  return-stroke  fields  (im/i/ni  AY-  l  ea  -  :  .  pr 

95S-958.  1 980 

[I7j  C  D  Vkeidmati  [  P  kodet.  and  M  A  l  man.  i  . 

amplitude  spectra  in  the  interval  trot:.  Hhi  kHz  n  2*1  S':*. 

iieophw  Res  let:  voi  s.  pp  9% I  4  -  -i .  lusi 
|  IK]  R  k  Baum.  '  Aunoine  lightning  characten/aiion  in  nchtti.iii. 
technoloev .  NASA  c'onf  Pub  2I2K.  I  A.A  Rl)  xu  30.  pp  i  l1J 
Apr  1980 

;|4j  1)  W  (1,'tord.  f  I'  kuder.  and  M  A  l  man.  A  ».im'  v! 

siitvmii.ii  sL\onil  rise  iniic  iigh'mng  curiem  [ui  ses  *.  ;  use  .i 
iiiduced  c oupltnc  siudies  presented  at  i*r:9  0  E|  l m  Svmp 
EMC  .  San  l>.cg<>  C  A  Oc?  1 479 

|  20 ]  E  L  Pitts  and  M  E  Thomas.  '  1 9S0  direct  strike  lightning  d;rtv. .  ” 
NASA  lech  M  emo  K|w4(*.  Langley  kesearc h  C  enter .  Hamp..  i 
VA-leb  |9S| 

1 2  1 1  \  L  Pitts  and  M  E  Thomas.  I9K I  Jtrcv  i  strike  liehtmne  data 

NASA  1  ech  Memo  k32~3.  Langley  Research  Center .  Hampton. 
\A  Mar  1982 

1 22 ]  M  A  l  man.  I)  k  McLain,  and  P.  P  kndgr.  The  eiecir. 

magnetic  radiation  trom  a  linue  antenna. "  Amer  J  Rh\±  .voi  4  V 
pp  33-38.  1975 

1 23]  Y  T  Lin.  M  A  Lnian.  and  R  B  Standler.  "Lightning  return 
stroke  models. "  d  Geophw  Res  .  v<d  85.  pp  |5"*|_|58V  Iv>Si' 

1 24]  M  J  Master.  M  A  Email  3  T  Lin.  and  R  B  stand. e? 

"Calculations  of  lightning  return  stroke  electric  and  magnetic  holds 
above  gnmnd.  J  iicttphw  Res  vol  Mfc.  pp  12127-12132.  J9SI 

|25]  M  A  l.'man.  l.tghtninK  New  Sork.  NY  McGraw  Hill.  1469 

[26]  M  A  l  man  and  E  P  Krtder.  "A  review  of  natural  lightning 

experimental  data  and  modeling.'1  / LEE  Trans  Llectroma$n 
Compat  .  vol  EMC  24.  pp  79-112.  1982 
1 27 1  k  Berger.  R  B  Anderson,  and  H  kronmger. •* Parameters «it 
lightning  flashes.**  Electro,  vol  41.  pp  23-3*7.  1975 
|28|  E  Garhagnati.  E  Guidice.  G  B  LoPipcro.  and  L  Magagnoh. 


40 


416 


IEEE  TRANSACTIONS  ON  ELECTROMAGNETIC  COMPATIBILITY.  VOL.  KMC  24.  NO  4.  NOVEMBER  196 


R 1 1 ic vi  dell  caratteristiche  dei  Futmmi  in  Italia  Risultati  ottcnuti 
negli  anni  1970-1973. “  I'Elettrotechmca.  vol  62.  pp  237-249 
1975 

{29J  A  J.  Eriksson.  “Lightning  and  tall  structures.  7> uns  South 
African  Inst  Elec  En%  .  vol  69.  pp  2-16.  Aug  I97H 
1 301  R  B.  Anderson  and  A.  J  Eriksson.  'Lightning  parameters  tor 
engineering  application. ’*  Electro .  vol.  69.  pp  65-102.  I9KU 


|3  1 J  F  Popolansk) .  'Frequency  distribution  ot  amplitudes  ot  iignming 
currents.  Electro,  vol  22.  pp  139-147.  Mj\  19*’? 

132}  D  W  Cliltord  and  H  W  Kasemir.  Triggered  lighimne  IEEE 
Trans  Electromunn  Compot  .  vol  KMC -24  pp  112-122  .'^2 
1*3)  K  H  Goldc.  Lightning  conductor.  in  LtKhtnine.  \  <>/'  II 
l.tchtniHc  Pnueitum  R  H  Cioldc.  td  New  York  \\  •Xc.idemic 
Press.  1977.  pp  545-576 


REFERENCES 


B.  Djebari,  J.  Hamelin,  and  C.  Leteinturier ,  Comparison  Between 

Experimental  Measurements  of  the  Electromagnetic  Field 
Emitted  by  Lightninq  and  Different  Theoretical  Models  - 
Influence  of  the  Upward  Velocity  of  the  Return  Stroke, 

Proc.  Electromagnetic  Compatibility  Conference,  Zurich, 
Switzerland,  March  1931. 

A.  J.  Ericksson,  The  CSIR  Lightning  Research  Mast  -  Data  for  1972- 
1982.  Internal  Report  No.  EK/0/82,  File  No.  E/73/2, 
National  Electrical  Engineering  Research  Institute,  Pre¬ 
toria,  South  Africa,  August  1982. 

R.  P.  Fieux,  C.  H.  Gary,  B.  P.  Hutzler,  A.  R.  Eybert-Bernard ,  P.  L. 

Hubers,  A.  C.  Meester,  P.  H.  Perroud,  J.  H.  Hamelin,  and 
J.  M.  Person,  Research  on  Artifically  Triggered  Liahtning 
in  France,  IEEE  Trans.  PAS,  PAS-97,  725-733,  1978. 

V.  I.  Idone,  and  R.  E.  Orville,  Liqhtning  Return  Stroke  Velocities 

in  the  Thunderstorm  Research  International  Program  (TRIP) , 
J.  Geophys.  Res.,  87,  4903-4915,  1982. 

D.  M.  Jordan,  and  M.  A.  Uman ,  Variation  in  Liaht  Intensity  with 

Height  and  Time  from  Subsequent  Lightning  Return  Strokes, 

J.  Geophys.  Res.,  88,  xxxx-xxxx,  1983. 

D.  M.  LeVine,  and  R.  Meneghini,  Simulation  of  Radiation  from  Light¬ 
ning  Return  Strokes:  The  Effects  of  Tortuosity,  Radio 
Science,  13,  801-809,  1978a. 

D.  M.  LeVine,  and  R.  Meneghini,  Electromagnetic  Fields  Radiated 

from  a  Lightning  Stroke:  Application  of  an  Exact  Solu¬ 
tion  to  Maxwell's  Equations,  J.  Geophys.  Res.,  83,  2377- 
2384,  1978b. 

M.  A.  Uman,  D.  K.  McLain,  R.  J.  Fisher,  and  E.  P.  Krider,  Currents 

in  Florida  Lightning  Return  Strokes,  J.  Geophys.  Res.,  78, 
3530-3537,  1973. 

C.  D.  Weidman,  and  E.  P.  Krider,  The  Fine  Structure  of  Liahtning 

Return  Stroke  Waveforms,  J.  Geophys.  Res.,  83,  6239-6247, 
1978  . 

C.  D.  Weidman,  and  E.  P.  Krider,  Submicrosecond  Risetimes  in  Light¬ 
ning  Return  Stroke  Fields,  Geophys.  Res.  Lett.,  7,  955-958, 
1980a. 

C.  D.  Weidman,  and  E.  P.  Krider,  Submicrosecond  Risetimes  in  Light¬ 
ning  Radiation  Fields,  in  Lightning  Technology,  NASA  Con¬ 
ference  Publication  2128,  FAA-RD-80-30 ,  1980b. 


42 


1 


C.  D.  Weidman,  E.  P.  Krider,  and  M.  A.  Uman ,  Lightning  Amplitude 

Spectra  in  the  Interval  from  100  kHz  to  20  MHz,  Geophys. 
Res.  Lett.,  8,  931-934,  1981. 

C.  D.  Weidman,  The  Submicrosecond  Structure  of  Lightning  Radiation 
Fields,  PhD  Dissertation,  University  of  Arizona,  1982. 


43 


<■  U.S.Govtrnment  Printing  Office  1984  759  062/644 


